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Abstract 

An  analytical  investigation  was  performed  to  evaluate 
the  influences  of  moisture  and  temperature  on  the 
bifurcation  load  of  cylindrical,  composite  panels  subject  to 
axial  loading.  The  composite  panels  ware  8-ply 
graphi te/epoxy  (AS/3501-5)  laminates.  The  laminate  ply 
orientations  considered  were  CO. , 45. , -45. , 90. 1^,  C90.,45., 
-45. ,0. 1_,  and  C45. , -45. , The  analysis  included  several 
different  panel  radii,  two  sets  of  boundary  conditions,  and 
three  sets  of  initial  moisture  conditions.  To  evaluate  the 
influences  of  moisture  and  temperature,  the  transverse 
moduli,  E»,  and  shear  moduli,  G,_.  were  degraded  based  on 
test  data  for  the  AS/3501-5  systF;m.  Each  ply  orientation 
was  evaluated  at  20  ti me/temperature  conditions  that  ranged 
from  80  F  to  300  F  and  moisture  concentrations  ranging  from 
zero  moisture  content  to  an  equilibrium  moisture 
concentration  distribution.  The  bifurcation  loads  were 
determined  using  the  STAGS-C1  finite  element  shell  analysis 
program.  The  bifurcation  analysis  mode  with  a  pre-buck  led 
linear  displacement  option  was  used  for  the  analysis. 
Moisture  and  temperature  were  found  to  cause  a  reduction  in 


the 

panels 

bifurcation 

load  rariqing  from 

the 

CO. ,45. , 

-45. ,90. 

laminate  to  42.7 

C45. 

• —45. 

laminate. 

MOISTURE  AND  TEMPERATURE  EFFECTS 


ON  THE 

INSTABILITY  OF  CYLINDRICAL  COMPOSITE  PANELS 


I.  Introduction 


Advanced  composites  are  a  leading  candidate  in  the 
search  -for  improved  aircraft  and  spacecraft  structural 
efficiency  and,  thus,  improved  system  performance.  In 
aircraft  applications,  until  very  recently,  the  use  of 
advanced  composites  has  beer,  limited  to  secondary  structural 
applications  or  unique  primary  structural  appi i cati ons. 
Howev  the  utilization  of  advanced  composites  for  aircraft 
primary  structure  has  begun  with  the  new,  all-composite, 
Lear  business  jet. 

The  desire  to  improve  the  efficiency  of  the 
structural  design  through  the  application  of  composites  in 
order  to  tailor  the  component's  strength  and  stiffness  to 
match  the  load  and  stiffness  requirements  results  in 
components  with  unique  structural  responses.  One  difficulty 
in  such  applications  of  advanced  composites  is  the  inability 
to  use  conventional,  classical  structural  design  practices 
to  predict  the  structural  response  of  composite  laminates. 
Attempts  to  apply  isotropic  analysis  techniques,  even  for 
quasi -isotropic  laminates,  yields  questionable  results. 
Often,  designers  and  stress  analysts  will  attempt  to 
compensate  with  the  use  of  special  knock-down  factors  to 


reduce  the  composite’s  strength  or  stiffness  so  that  a 


“conservative"  design  is  produced.  The  efficient  use  of 
advanced  composites  in  primary  structural  applications  will 
depend  upon  the  development  of  suitable  design  and  analysis 
tools  that  adequately  predict  the  composite’s  structural 
response. 

In  a  semi-monocoque  structural  design,  thin  skin  panels 
along  with  fuselage  frames  and  longerons  or  wing  skins  with 
spars  and  ribs  are  used  to  achieve  an  efficient  structural 
design.  Several  al 1 -composite  aircraft  designs  have  been 
proposed  which  utilize  a  semi-monocoque  design  approach  in 
which  composite  panels  are  used  for  the  fuselage  skins  or 
wing  skins  and  composite  members  are  used  for  the  frames, 
longerons,  etc. 

In  such  a  semi-monocoque  design,  the  accurate 
determination  of  the  buckling  load  of  the  curved  skin  panels 
is  necessary  to  prevent  premature  structural  failure.  This 
determination  of  the  buckling  loads,  even  for  relatively 
simple  isotropic  structures,  is  often  complicated  by  unusual 
shapes  and  boundary  conditions  and/or  combinations  of  loads 
and  moments.  Designers  and  analysts  rely  upon  simple 
theoretical  cases  with  known  solutions  or  parametric 
equations  based  upon  laboratory  testing  to  develop  buckling 
allowables.  The  application  of  numerical  techniques, 
especially  finite  difference  and  finite  element  techniques, 
has  provided  a  new  method  of  performing  shell  stability 
analyses. 
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These  numerical  techniques  are  of  special  importance  in 
the  evaluation  of  composite  curved  panels.  Unlike  isotropic 
panel 5,  the  uniqueness  of  the  composite  panel,  due  to  the 
design  dependent  number  and  orientation  of  the  plies, 
prevents  the  development  of  general  design  tools.  Hence,  in 
most  applications,  these  numerical  techniques,  in  the  form 
of  computer  programs,  provide  the  only  means,  besides 
experimentation,  for  determining  a  composite  panel's 
buckling  load. 

Perhaps  because  of  the  current  limited  applications  of 
composites  for  thin,  structural  skin  panels,  there  has  been 
very  little  published  on  the  stability  of  curved  composite 
panels.  The  majority  of  work  has  dealt  with  flat  plates  or 
shel ls-of-revolution.  J.  D.  Wilkins,  Ref  1,  conducted  an 
experimental  investigation  of  the  buckling  strength  of  a 
cylindrical  composite  panel  subject  to  axial  compression 
with  two  sets  of  boundary  conditions.  M.  L.  Becker  (2) 
extended  Wilkins'  work  to  more  boundary  conditions,  several 
panel  sizes  and  aspect  ratios,  and  different  ply 
orientations.  Becker  also  used  the  finite  difference 
version  of  the  STAGS-C  (Structural  Analysis  of  General 
Shells)  program  to  compare  the  analytical  and  experimental 
buckling  results. 

An  aspect  of  composite's  research  receiving 
considerable  attention,  is  the  influence  of  moisture  and 
temperature  on  the  structural  characteri sties  of  the 
composite  system.  These  influences,  sometimes  referred  to 
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as  hygrothermal  effects,  fall  in  the  general  category  of 
environmental  influences.  Such  environmental  effects  have 
been  found  to  significantly  degrade  the  mechanical 
properties  of  most  organic  matrix  compounds  <3-14). 

Purpose 

This  thesis  extends  the  analytic  work  of  Becker  to 
evaluate  the  instability  of  composite  cylindrical  panels 
subject  to  axial  compression,  using  a  new  finite  element 
version  of  the  STAGS-C  shell  analysis  program.  The 
influence  of  environmentally-caused  reductions  in  the 
composite’s  transverse  and  shear  moduli  on  the  panel’s 
bifurcation  load  will  be  investigated  using  this  finite 
element  program. 

Scope 

The  cylindrical  panels  evaluated  in  this  thesis  are  12 
in.  by  12  in.,  3-piy  laminates  made  of  Graphite/Epoxy.  The 
panel's  thickness,  width,  and  height  are  held  constant  while 
the  radius,  ply  orientations,  and  boundary  conditions  are 
varied.  Three  radii,  three  ply  orientations,  and  two 
boundary  conditions  are  investigated.  During  the 
investigation  of  the  environmental  influences,  four 
temperatures  and  three  surface  moisture  conditions  at  five 
time  values  are  evaluated. 


IL  Environmental,  Influences  on  Composites 


Moisture  and  Temperature  Influences 

While  composites  have  many  superior  properties  compared 
to  metals,  the  commonly  used  epoxy  composite  systems  are 
significantly  affected  by  environmental  exposure  to  an  agent 
such  as  water,  which  is  absorbed  by  the  polymer  resin,  and 
to  thermal  conditions  which  are  near  or  exceed  the  polymer’s 
glass  transition  temperature.  The  fibers,  which  are 
typically  graphite  or  boron,  are  not  affected  by  either 
water  or  the  moderate  temperatures  an  aircraft  would 
encounter  during  normal  service.  As  a  consequence,  the 
resin-dominated  material  properties  are  affected. 

Whitney  and  Ashton  (3)  reported  that  the  environmental 
factors  which  influence  the  resin  properties  were  a) 
increases  in  temperature,  b)  absorption  by  the  polymeric 
resin  material  of  a  swelling  agent  such  as  water  vapor,  and 
c)  the  sudden  expansion  of  absorbed  gases  in  the  resin.  It 
was  noted  that  resin  swelling  alone,  due  to  moisture  and 
temperature  effects,  could  cause  a  flat  composite  plate  to 
buckle. 

Temperature  and  moisture  act  on  the  polymeric  resin  in 
two  ways.  Increasing  temperatures  cause  the  resin  to  swell. 
Absorbed  water  vapor,  which  in  one  study  (4)  was  believed  to 
be  due  to  the  water  molecules  oonding  to  the  hydroxyl  groups 
in  the  epoxy  polymers.  also  causes  the  resin  to  swell. 
Matrix  swelling  and  rapid  heating  may  eventually  lead  to 
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surface  crazing  and  surface  cracks  which  will  affect  the 
resin’s  mechanical  properties  (12  and  13). 

Temperature  and  moisture  also  act  together  to  cause 
reductions  in  the  temperature  range  over  which  resin 
mechanical  properties  are  fairly  stable.  Moisture  absorbed 
in  the  resin  results  in  a  plasticization  of  the  resin.  This 
plasticization  is  the  result  of  the  lowering  of  the  glass 
transition  temperature  T  .  The  glass  transition  temperature 
is  actually  a  temperature  range  in  which  below  this  range 
the  resin  is  essentially  brittle  and  above  this  range  it 
behaves  rubbery. 

This  change  in  T  ,  due  to  an  absorbed  agent,  is  common 
in  polymer  science.  The  polymer’s  Tg  is  a  function  of  the 
free  volume  in  the  material.  It  is  accepted  theory  (14) 
that  1/40  of  the  total  volume  of  the  material  is  free  volume 
(total  volume  less  the  molecular  volume)  at  and  below  the 
material’s  T  .  As  the  polymer  absorbs  an  agent,  such  as 
water,  which  contains  more  free  volume,  the  percentage  of 
free  volume  of  the  mixture  will  increase  (12).  When  the 
mixture’s  free  volume  is  above  1/40  of  the  total  volume, 
the  resin  begins  to  soften.  Thus,  the  mixture’s  temperature 
must  be  lowered  for  the  mixture’s  free  volume  to  return  to 
1/40  of  the  total  volume.  This  results  in  a  lower  T  .  As 
increasing  concentrations  of  moisture  are  absorbed  in  the 

resin,  the  T  is  continual lv  lowered  until  a  moisture 

y 

equilibrium  is  reached.  Once  equilibrium  is  reached,  the  T^ 
remains  constant. 
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These  changes  in  the  resin  have  been  found  to  result  in 
decreases  in  the  tensile  properties  (6)  and  reductions  in 


the  transverse  and  shear  moduli  (7,  8,  and  13)  of  the 
composite.  A  slight  increase  in  the  longitudinal  elastic 
moduli  was  reported  in  Ref  13.  As  noted  in  Ref  12,  these 
changes  in  the  composite’s  mechanical  properties  may  be 
grouped  into  two  general  classifications:  (a)  those  changes 
due  to  moisture-induced  plasticization  which  reduces  the 
temperature  range  over  which  the  material  properties  are 
relatively  stable  and  <b)  those  losses  due  to  mechanical 
damage  from  moi sure-i nduced  swelling  and  rapid  heating  which 
results  in  surface  crazing  and  cracking  which  affect  the 


mechani cal 
mechanical 


properti es 
properties 


below  the  T  .  The  reduction  in 
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is  most  significant  in  the  tensile 


properties.  This  is  primarily  because  the  resin  is  very 


brittle  below  the 
shear  to  any  surface 


T  and  the  resin  is 

g 

flaws  than  it  is  in 


less  sensitive  in 
tension  ( 13)  . 


of  Absorbed  Moisture 

When  the  mechanical  properties  are  determined  by 
testing,  they  are  usually  measured  at  known  temperatures  and 
moisture  concentrations.  The  absorption  of  water  vapor  into 
the  resin  is  through  diffusion.  Under  certain 
circumstances,  Pick’s  second  law  of  diffusion  (23)  has  shown 
good  correlation  with  test  data.  Fig  1,  at  predicting  the 
rate  of  moisture  weight  gain  in  a  polymeric  resin  and  in  the 
composite.  Pick  developed  this  equation  in  1855  by  drawing 


Figure  1.  Moisture  Gain  as  a  Function  of  Time  for  AS/3501-5  Composites 
(Ref  13) 


l 
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an  analogy  between  heat  conduction  in  a  solid  and  diffusion 
through  a  solid. 


The  Fick  equation  isi 

ac  _  y  a2c 
at  az2 

where!  C  -  the?  concentration  of  the  diffusing 
substance-  through  the  thickness  of 
the  laminate  as  a  function  of  time 
and  distance  through  the  thickness 

Z  -  the  space  coordinate  measured  normal 
to  the  surface 

K  -  the  diffusion  constant 

t  -  time 

The  solution  of  this  partial  differential  equation  with 
boundary  and  initial  conditions  pertinent  to  the  problem  is 
shown  below.  This  series  solution  in  a  slightly  different 
form  is  found  in  Section  4.3.3  of  Pef  23. 


CU,t)~Cx+(Ci-C.)Z 

h 

CX>  ,1/ 

.  Z  V*  C*  cos  nTC  -  Cx  Sin  ivscZ  c%p  —  rz — 
TC  L  n  h  K 

H.iLfiaf’  exp 

1C  i-i  zm-ri  h  h*" 

tn«i 
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C0,  Cj,  and  represent  the  inital  set  of  moisture 
conditions  for  the  laminate.  is  the  initial  moisture 
concentr ati on  distribution  through  the  thickness  of  the 
laminate.  Cj  is  the  initial  moisture  condition  at  the 
surface  of  one  side  of  the  laminate.  C0  i s  the  initial 
moisture  condition  at  the  surface  on  the  other  side  of  the 
laminate.  The  thickness  of  the  laminate  is  represented  by 
h. 

Using  this  series  solution  with  a  known  diffusion 
constant  and  prescribed  initial  conditions,  the  moisture 
concentration  distribution  through  the  thickness  can  be 
determined.  With  the  assumption  that  the  effective  moisture 
concentration  of  each  ply  can  be  approximated  by  the 
calculated  moisture  concentration  at  the  middle  of  the  ply, 
the  reduced  mechanical  properties  of  each  ply  can  be 
determined  from  appropriate  test  data. 

This  series  solution  is  seen  to  represent  the  summation 
of  a  steady  state  moisture  distribution,  represented  by  the 
first  and  second  terms,  and  a  transient  moisture 
distribution  that  is  a  function  of  time.  This  is 
represented  by  the  last  two  terms.  The  influence  of  the 
transient  terms  decreases  with  increasing  time.  The 
accuracy  of  this  series  approximation  is  dependent  upon  the 
number  of  terms  used  during  the  two  summations.  This 
solution  was  calculated  by  the  computer  program  shown  in 
Appendix  A.  To  insure  an  accurate  solution  the  summations 
were  carried  out  until  there  was  no  change  from  the  previous 
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answer 


On  the  CDC  computer  this  equated  to  14  significant 


digits  of  accuracy. 

There  are  limitations  on  the  application  of  Fick’s 
equation.  The  series  solution  of  Fick’s  equation  assumes 
that  the  moisture  diffusion  coefficient  K.  is  constant.  The 
diffusion  coefficient  is  a  function  of  the  laminate’s 
temperature.  However,  since  the  moisture  diffusion  is  a 
relatively  slow  process,  with  many  months  or  years  required 
before  the  moisture  ccnccntrati on  distribution  through  the 
laminate  achieves  equilibrium!,  the  diffusion  process,  in 
simple  cases,  may  be  assumed  to  take  place  at  a  constant 
temperature.  Eergmann  and  Nitsch  (4)  have  noted  that  K.  also 
varies  with  the  laminate’s  moisture  concentration  generally 
increasing  with  increasing  moisture  concentrat i on  levels. 

The  accuracy  of  Fick’s  equation  is  also  affected  by 

rapid  temperature  changes.  Rapid  thermal  heating  of  the 

laminate  (which. may  be  due  to  flight  at  supersonic  speeds), 

where  the  laminate  is  heated  to  temperatures  near  the 

material’s  T  ,  has  been  found  to  increase  the  rate  of 
9 

moisture  weight  gain  above  that  predicted  by  the  Fick 
equation  (7,  10,  12,  and  13).  This  increase  is  believed  due 

to  the  development  of  surface  crazing  and  cracking  brought 
about  by  the  rapid  heating  and  resin  swelling  (13). 

With  the  restrictions  of  no  rapid  heating  and  no 
surface  crazing  or  cracking,  and  assuming  that  K  is 
constant,  Fick's  equation  has  been  generally  accepted  as  a 
good  initial  approximation  of  the  moisture  concentration 
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distribution  for  simple  cases  <4,  7,  9,  11,  and  13). 


AS/3501-5  h-^chan  i  cal^  Properties 

The  GTAGS-C1  shell  analysis  program  discussed  in 
Section  III  requires,  as  input  parameters,  the  composite’s 
longitudinal  moduli  E^,  transverse  moduli  E^,  shear  moduli 
G  and  Poisson’s  ratio  V  .  Poisson’s  ratio  V  relates 
the  strain  in  the  2  direction  to  the  strain  in  the  1 
direction  when  stresed  in  the  1  direction.  Experi mental 1 y- 
measured  data  for  a  graphite/epoxy  system,  AS/3501-5,  from 
Fig  S. 18  of  Ref  15  was  used  in  the  determination  of  the 
elastic  moduli  as  a  function  of  temperature  and  moisture 
concentration.  The  units  of  stress  in  GPa  and  temperature 
in  degrees  Kelvin  are  converted  to  psi  and  Fahrenheit.  The 
values  of  E_  and  G^  used  in  this  work,  from  which 
intermediate  values  were  interpolated,  are  shown  in  Table  I. 
The  method  used  for  i ntorpol ati on  is  described  in  the  next 
section.  Figure  2  shows  this  data  in  a  graphical  form. 

The  moisture  and  temperature  influences  on  the 
transverse  and  shear  moduli  are  clearly  evident  in  the 
experimental  data  for  AS/3501-5  shown  in  Table  I  and  Fig  2. 
The  transverse  moduli  E^  shows  deg-adati on  both  at  room 
temperature  and  at  elevated  temperatures  while  the  shear 
moduli  G^  only  shows  degradation  at  elevated  temperatures. 
The  moisture  caused  chanus  in  the  T  and  the  resulting 
plasticization  of  the  resin  is  shown  by  the  increased 
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degradation  in  the  moduli  with  increasing  moisture 
concentration  at  each  elevated  temperature.  This  is 
consistent  with  the  previously  noted  expected  changes  in 
elastic  moduli.  The  longitudinal  moduli  is  dominated  by 
the  -fiber  stiffness  and  hence  is  not  significantly  influeced 
by  changes  in  moisture  and  temperature  as  are  the  matrix 
dominated  E^  and  moduli. 

The  value  of  used  in  this  analysis  was  taken  from 
Table  1.9  of  Ref  15.  The  value  of  V"12  0.30  which  is  a 
representat i ve  value  for  AS/3501-5  composites.  Examination 
of  the  test  data  (13)  upon  which  the  material  in  Ref  15  was 
based  shows  that  V"12  ^oes  vary  with  temperature  and 
moisture  concentration.  The  range  of  variation  was  from 
0.36  to  0.46.  The  sensitivity  of  the  analyses  to  this 
variation  in  "l^2  was  found  to  be  very  small.  This  is 
discussed  in  more  detail  in  Section  IV. 

The  experimental  stress-strain  data  from  Ref  13, 
reflecting  the  changes  in  the  shear  and  transverse  moduli  as 
a  function  of  temperature  and  moisture,  are  shown  in  Figs  3 


and  4,  respectively 


Table  I 


Values  of  Transverse  and  Shear  Moduli 
for  AS/3501-5  (15) 

Transverse  moduli,  E_,  <psi) 


Moisture  Temperature 

Concentr  at i on 


(percent ) 

SO  F 

200  F 

250  F 

300  F 

0.0 

1 . 41375E06 

1.09473E06 

1.015E06 

1.015E06 

0.050 

1 . 305E06 

0.9135E06 

0. 6235E06 

0.522E06 

1.050 

1.2615E06 

0.B41E06 

O.4785E06 

0.290E06 

Shear  Modul i , 

Gl2  <psi> 

0.0 

0. 8555E06 

0. 7830E06 

0.6S15E06 

0. 6525E06 

1.050 


0 . 8555E06  0 . 6597E06 


0.391 5E06  0 ..  1 5225E06 


Wet  =  Equilibrium,  7 55?  R.H.  @  160  F 


- , - t - , - , - , 

0  1.0  2.0  3.0  4.0  5.0 

Shear  Strain  -  7l2  (Percent) 

Figure  3«  Typical  In- plane  Shear  Stress-strain  Curves  for  AS/3501-5 
Unidirectional  Composites  (Ref  13) 


1  0  0.2  0.4  0.6  0.6  1.0 

'  Strain  (Percent) 

i 

i  Figure  4.  Typical  Transverse  Stress-strain  Curves  as  a  Function  of 

Temperature  for  A5/3501-5  unidirectional  Composites  (Ref  13) 
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Hit  Analytical  Approach 


Background 

The  cylindrical  shell  was  one  of  the  first  shell  shapes 
to  be  used  structurally.  From  an  analytical  viewpoint,  the 
analysis  of  cylindrical  shells  and,  in  particular,  the 
analysis  of  the  stability  of  cylindrical  shells  was  begun  in 


the  late 

1800s 

.  In  the  early 

1900s,  the 

analytical 

treatment 

of 

cylinders  under  a 

variety  of 

loads  was 

developed 

with 

the  most  comprehens 

ive  treatment 

done  by 

Flugge  in  1932.  In  1933,  Donnell  (16)  proposed  a  set  of 
simplifications  of  the  non-linear  cylindrical  shell 
equilibrium  equations.  These  Donnell  equations  formed  the 
basis  for  much  of  published  stability  analyses  until,  in 
recent  years,  the  availability  of  high  speed  computers  has 
led  to  the  development  of  efficient  numerical  analytical 
techniques. 

The  application  of  numerical  techniques  to  structural 
problems  began  in  the  1950s  with  the  advent  of  matrix 
analysis.  These  matrix  techniques  have  evolved  through  the 
develooment  of  a  strong  mathematical  base  into  the  familiar 
finite  difference  and  finite  element  methods  now  widely 
used.  In  the  application  of  these  techniques  to  the  study 
of  the  stability  of  shells,  a  series  of  programs  were 
developed  at  the  Lockheed  Palo  Alto  Research  Laboratory. 
These  programs  have  generally  been  known  bv  the  name  STAGS. 
The  two  most  recent  versions  of  the  program  are  the  finite 
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difference  version,  STAGS-C,  and  the  finite  element  version, 


STAGS--C1.  M.  Becker  used  the  finite  difference  version, 
STAGS-C,  in  his  comparison  of  experimental  and  analytical 
buckling  loads  for  cylindrical  composite  panels  (2).  Since 
his  work,  the  finite  element  version,  STAGS-C1,  has  been 
released  and  was  used  in  this  thesis. 

Aeelicatign  t?  §13511  Analysis 

The  STAGS-C1  program  has  been  specifically  developed  to 
perform  a  collapse  analysis  of  general  stiffened  and 
unstiffened  shells.  The  program  has  several  operating 
modes.  Among  these  are  a  linear  or  geometric  non-linear 
static  analysis,  a  bifurcation  analysis  with  a  linear  or 
non-linear  stress  state,  and  a  small  vibration  analysis 
with  either  a  stress-free,  linear,  or  nonlinear  stress 
state.  The  program  is  capable  o*  handling  both  isotropic  or 
layered  orthotropic  materials.  The  reader  is  referred  to 
Ref  17  for  a  comprehensive  review  of  the  capabilities  of 
this  program  and  to  Ref  18  for  an  explanation  of  the 
program’s  input  parameters. 

Both  the  finite  difference  and  finite  element  programs 
are  useful  in  the  analysis  of  shells.  While  the  STAGS-C 
finite  difference  program  is  usually  considered  better 
suited  for  modeling  curved  surfaces  due  to  the  curvilinear 
fit,  it  has  been  found  to  be  less  efficient,  in  certain 
analysis  mooes,  than  the  STAGS-C1  finite  element  program 
which  uses  flat  elements.  Curved  finite  elements  have  not 


yet  been  found  to  be  practical  because  they  tend  to  require 
large  amounts  of  computer  time  for  the  formulation  of  the 
stiffness  matrix.  They  also  have  the  problem  of  self¬ 
straining;  straining  do  to  rigid  body  movement  associated 
with  the  out-of-plane  displacement  of  the  shell’s  surface 
prior  to  buckling  (19).  Until  suitable  curved  elements  are 
developed,  modified  flat  elements  are  used  in  STAGS-C1. 

The  modeling  of  a  cylindrical  surface  with  flat 
elements  is  shown  in  Fig  5.  These  flat  elements  have 
conformity  problems  which  are  of  importance  in  their 
application  to  nonlinear  and  stability  problems.  As  seen  in 
Fig  5,  both  rotational  and  displacement  problems  develop  at 
the  nodes.  The  rotational  problem  is  dealt  with  by  assuming 
that  the  angle  of  intersection  <*■  is  small,  and,  as  a 
consequence,  the  normal  rotation  is  ignored  and  the 
conformity  constraint  is  met  by  letting/9^1*  =/?  . 

The  interelement  displacement  continuity  is  more 
difficult  to  deal  with.  Complete  displacement  compatabi 1 i ty 
along  the  common  boundary  requires  that 

-  v lz) )  COS  ( */2)  -  (w (,)+  TV (z))  sm  ( Vz)- o 
(Y/(t)-WiZ))  COS  (°t/2>  +  (ir(,Kv{Z))  SIR  (*/z)  -  o 

where  v  and  w  are  the  displacements  in  the  y  and  z 
directions,  respectively,  as  shown  in  Fig  5.  If  v  and  w  are 
not  represented  by  polynomials  of  the  same?  order,  then  these 
conditions  will  not  be  met.  Failure  to  meet  these 
conditions  will  allow  the  individual  elements  to  buckle 
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under  axial  compression  as  plates  with  -free  edges 


This  displacement  conformity  is  met  by  adding 

additional  freedoms  to  the  element  that  will  raise  the  order 

of  the  polynomials  representing  inplane  deformation.  Since 

w  is  represented  by  a  cubic  polynomial,  because  the  strain 

energy  expression  includes  second  order  derivatives  of  the 

transverse  displacement  w,  it  is  necessary  that  u  and  v  also 

be  represented  by  cubic  polynomials.  This  is  achieved  by 

the  use  of  two  rotations  at  each  node,  -v,  and  u,  ,  and 

x  y 

tangential  displacements  at  midside  nodes.  The  difference 
between  these  two  rotations  is  the  shear  strain  at  eacn 
node.  This  shear  strain  ^2’  ^3’  ^4*  *s  in*rcduced  as 

an  additional  freedom  .  Thus,  each  element  has  7  degrees- 
of-freedom  at  each  node  and  4  additional  displacement 
freedoms  for  a  total  of  32.  These  additional  degrees-of- 
freedom  were  added  to  raise  the  order  of  the  polynomial 
representing  the  inplane  deformations.  This  element  is 
referred  to  in  the  STAGS-C1  program  as  SH411. 

A  somewhat  simpler  version  of  the  SH411  element,  the 
SH410  element,  is  also  included  for  thin  shell  analysis. 
This  element  excludes  the  midside  tangential  displacements 
and  uses  only  an  average  normal  rotation  at  each  cf  the 
corner  nodes.  This  restricts  u  to  a  linear  function  in  the 
x-direction  and  v  to  a  linear  function  in  the  y-direction. 
The  shear  strain  is  also  suppressed  at  the  nodes. 

The  freedoms  for  both  these  elements  are  shown  in  Fig 
6.  The  reader  is  referred  to  Refs  19  and  20  for  a  more 
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thorough  explanation  of  the  development  of  these  elements 
Reference  19  contains  the  derivation  of  the  SH411  element. 


Classical  Laminated  Plate  Theory 

Since  this  thesis  deals  with  the  problem  of  the 
buckling  of  thin,  composite  panels,  a  statement  of  the 
pertinent  assumptions  and  restrictions  is  needed.  Reference 
21  contains  a  concise  statement  of  restrictions  and 

assumptions  applicable  to  thin  plate  theory. 

Restrictions! 

A)  Each  layer  of  the  laminate  is  orthctropic,  linear 
elastic,  and  of  constant  thickness. 

B)  The  panel  thickness  is  very  small  compared  to  its  length 
and  width. 

C>  The  radius  of  curvature  of  the  panel  is  large  insuring 

that  the  curved  panels  may  be  described  as  shallow  shells  or 

quasi -shal 1 ow. 

Assumptions: 

A)  Stresses  acting  in  the  xy  plane  dominate  the  panel 
behavior.  As  a  consequence,  <TZ ,  Txz,  and  Tyz  are  assumed  to 
be  zero  and  a  state  of  plane  stress  exists. 

8)  The  Kirchhoff  assumption  of  negligible  transverse  shear 
strains,  $  xz  ^yz’  anc*  transverse  normal 

strain,  applies. 

A. 

C)  The  displacements  u,  v,  and  w  are  small  compared  to  the 
panel  thickness. 

D)  The  square  of  the  rotations  are  small  compared  with  the 
linear  strains, 

E)  The  strains,  €  ,  £  ,  and  ^  ,  are  small  compared  to 

x  y '  x  >' 

unity. 
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A  diagram  of  the  plate  element’s  geometry  and  stress 
resultants  is  shown  in  Fig  7. 

The  assumption  that  the  transverse  shear  strains  may 
be  assumed  to  be  zero  results  in  the  transverse  shear 
stresses  also  being  zero.  This  is  a  standard  assumption  in 
classical  laminated  plate  theory.  This  is  fcrtunate  because 
the  SH411  and  SII410  plate  elements  described  previously  do 
not  include  transverse  shear  strains  in  their  formulation. 
However,  in  addressing  the  problem  of  environmental 
degradation  in  composites,  this  assumption  may  not  be  valid. 
Flaggs  and  Vinson  (22)  have  developed  a  general  buckling 
theory  for  flat  plates  which  accounts  for  the  moisture  and 
temperature  effects  and  includes  transverse  shear,  normal 
deformation,  and  bendi ng-oxtensi onal  coupling.  The 
application  of  this  theory  to  curved  panels  will  require  a 
finite  element  with  transverse  shear  capability.  One  such 
element,  referred  to  as  an  AHMAD-type  element  (19),  is 
currently  under  development  but  is  not  yet  available  in 
STAGS— Cl . 

Because  of  the  apparent  need  for  a  transverse  shear- 
capable  element  to  totally  address  the  moisture  anu 
temperature  effects,  the  study  of  these  effects  will  be 
limited  to  the  degredati ons  in  elastic  m.odul  i  due  to 
moisture  and  temperature.  Resin  swelling,  either  due  to 
moisture  absorption  or  temperature  expansion,  will  not  be 
included. 

Another  restriction  needed  to  use  STAGS-C1  is  that  each 
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layer's  material  properties  are  assumed  to  be  constant 
across  the  layer’s  thickness  and  throughout  that  layer.  In 
the  calculation  of  the  degraded  elastic  moduli,  each  layer 
was  assumed  to  have  a  uniform  moisture  concentrati on  equal 
to  the  value  calculated  at  the  ply’s  mid-surface.  In  the 
determination  of  the  temperature  degradation,  the  entire 
laminate  was  assumed  to  be  at  the  specified  temperature. 

With  these  additional  restrictions  and  assumptions,  the 
STAGS— Cl  program  should  adequately  predict  the  panel's 
critical  load.  By  limiting  the  problem  to  only  changes  in 
elastic  moduli,  the  influence  of  these  moduli  degradations 
may  be  addressed  separately.  Also,  the  problem  of  modeling 
the  laminated  plate  is  simplified. 

Because  flat  elements  are  being  used  in  the  finite 
element  model,  the  formulation  of  the  stress-strain 
relationships  is  straightf orward  as  shown  in  Ref  21.  Since 
unsymmetric  moisture  concentration  distributions  through  the 
laminate’s  thickness  will  be  investigated,  bending- 
extensional  coupling  will  develop.  The  equations  relating 
the  stress  resultants  to  the  mid-surface  strains  and 


curvatures  are 


where!  N.  —  membrane  stress  resultants 
1 

-  bending  stress  resultants 
Aij’  Bij’  Dij  “  stiffnesses 

€?  -  mid-surface  extensional  strains 
v  O 

0  -  mid-surface  shear  strains 

w  xy 

-  mid-surface  curvatures 

The  A.  .  are  called  the  extensional  stiffnesses,  the  B.  .  are 
i  J  i  J 

called  the  coupling  stiffnesses,  and  the  are  called  the 

bending  stiffnesses.  The  j  results  in  coupling  between 
bending  and  extension.  For  symmetric  laminates,  the  ^ 
vanish  and  there  is  no  coupling  between  bending  and 
extension.  The  influence  of  the  B^ ^  will  bo  evaluated 
through  the  appropriate  selection  of  initial  moisture 
conditions  that  result  in  symmetric  and  unsymmetric 


1 aminates 


Finite  Element  Model 


This  thesis  will  evaluate  the  stability  of  cylindrical 
composite  panels  subject  to  axial  compression.  A  typical 
composite  fuselage  skin  panel  with  backup  structure  is  shown 
in  Fig  8.  This  structural  component  is  used  as  the  basis 
for  the  development  of  the  finite  element  model  shown  in  Fig 
10.  The  STAGS-C1  cylindrical  shell  geometry  is  shown  in 
Fig  9. 

The  panel  is  square  with  the  panel  height  and 
circumferential  width  being  equal.  Such  a  panel  is  defined 
to  have  an  aspect  ratio  of  1.  The  aspect  ratio  is  defined 
to  be  equal  to  the  height  divided  by  the  width  for  both 
panels  and  elements.  All  panels  in  this  thesis  have  aspect 
ratios  of  1.  This  panel  aspect  ratio  determination  is 
different  from  Becker’s  (2)  where  the  panel’s  chord 
dimension  was  defined  as  the  panel’s  width.  This  resulted 
in  the  panel’s  aspect  ratio  being  a  function  of  the  panel’s 
radius  of  curvature.  The  use  of  this  notation  would  have 
made  the  comparison  of  panels  with  different  radii  difficult 
because  each  panel  would  have  a  different  aspect  ratio. 

The  boundary  conditions  used  in  the  evaluation  of  the 
moisture  and  temperature  effects  were  selected  to  represent 
those  of  a  typical  backup  structure  rather  than  the  standard 
simple  or  fixed  boundary  conditions.  Two  sets  of  boundary 
conditions  were  chosen.  The  first  set  of  conditions  assumes 
that  the  backup  structure’s  ring  frames  and  longerons  ace 


effective  in  restraining  out-of-plane  deflections  <W.'  and 
rotational  movement  (RU,  RV,  RW) .  The  top  frame  and  side 
longerons  are  assumed  to  not  be  effective  in  resisting  the 
axial  compressive  load  but  do  resist  in-plane  displacements 
in  the  circumf erenti al  direction  (V> .  The  second  set  of 
boundary  conditions  is  identical  to  the  first  with  the 
exception  that  the  longerons  are  assumed  to  have  no 
torsional  stiffness  and,  hence,  cannot  resist  a  torsional 
rotation.  In  the  first  set,  the  longerons  were  assumed  to 
be  closed  sections  and,  hence,  effective  in  torsion.  In  the 
second  set,  the  longerons  were  assumed  to  be  open  sections 
which  have  little  effective  torsional  resistance.  The 
boundary  conditions  are  summerized  in  Table  II. 

Table  II 

Panel  Boundary  Conditions 


U 

V 

w 

RU 

RV 

RW 

U 

V 

w 

RU 

RV 

RW 

Top 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

Right  Side 

1 

0 

0 

0 

0 

0 

1 

0 

0 

1 

0 

0 

Bottom 

0 

0 

c 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Left  Side 

1 

0 

0 

0 

0 

0 

1 

0 

0 

1 

0 

0 

where  0  represents  a  fixed  displacement  and  1  represents  a 
free  displacement  along  the  panel’s  edge. 

The  loads  applied  to  the  panel  represent  axial 
compression  leading.  This  load  is  assumed  to  act  uniformly 
across  the  top  of  the  pant?]  and  is  reacted  across  the  bottom 
of  the  panel . 
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Bifurcation  Analysis  Method 


STAGS-C1  has  two  buckling  analysis  modes.  One  uses  a 
linear,  pre-buckling  displacement  state  and  is  referred  to 
in  the  STAGS-C1  manual  as  the  bifurcation  analysis  with  a 
linear  stress  state.  The  second  method  uses  a  geometric 
pre-buckl ing,  nonlinear  displacement  calculation  and  is 
referred  to  as  the  bifurcation  analysis  with  a  non-linear 
displacement  state.  As  noted  in  Ref  19,  the  SH411  and  SH410 
elements  have  not  produced  good  results  using  the  nonlinear 
mode  in  predicting  post-buckling  modes.  Thus,  the 
bifurcation  analysis  with  a  linear  displacement  state  prior 
to  bifurcation  will  be  used.  This  method  will  calculate  the 
pre-buckling  displacements  and  rotations,  stress  resultants, 
strains,  and  stresses  as  desired.  It  also  predicts  the 
bifurcation  eigenvalue  and  the  shape  of  the  eigenvector.  It 
does  not  yield  any  post-buckling  information. 
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§i¥*iy£i£i°Q  of  STAGS-Cl  Program 


A  series  of  comparison  runs  with  the  STAGS-C  finite 
difference  version  of  the  program  and  a  convergence  study 
to  determine  the  best  mesh  size  was  performed.  Previously 
run  STAGS-C  results  for  a  12  in.  high  by  12.567  in.  wide 
cylindrical  composite  panel  (panel  aspect  ratio  cf  1  per 
Becker’s  notation),  with  several  sets  of  boundary 
conditions,  were  compared  with  the  equivalent  STAGS-Cl 
finite  element  results.  These  STAGS-C  finite  difference  runs 
were  done  by  Backer  as  a  part  of  the  work  reported  in  Ref 
4.  The  STAGS-C  runs  used  the  linear  bifurcation  analysis 
and  the  STAGS-Cl  runs  used  the  bifurcation  analysis  with  the 
linear  stress  state  option. 

The  panel's  ply  orientation  was  C45.,-45. with  a 
panel  radius  of  12  in.  The  boundary  conditions  along  the 
vertical,  straight  side  were  changed  while  the  boundary 
conditions  along  the  top  and  bottom  were  not  changed.  The 
STAGS-Cl  results  compared  well  with  the  STAGS-C  results  even 
though  the  STAGS-C  program  was  using  a  curved  surface 
formulation  and  the  STAGS-Cl  program  was  using  flat  plate 
elements.  The  results  of  this  comparison  study  are  shown  in 
Table  III.  These  STAGS-Cl  runs  were  made  using  the  simpler 
SH410  flat  plate  element  with  a  24  row  by  25  column  element 
mesh. 
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Table  III 


Comparison  of  STAGS-C  and  STAGS-C1  Results 


Bound. 

Case 

STAGS-C 

STAGS-C 1 

Percent 

Cond. 

Bif.  Load 

Bif.  Load 

Diff . 

Free 

BF1 

116.5 

107.2 

8.0 

CC1 

BC1 

392.  7 

399.0 

1.6 

CC4 

BD1 

498.0 

488.5 

1.9 

SSI 

BS1 

391.6 

396.  1 

1.1 

SS4 

BT1 

477.2 

452.6 

5.2 

When  using  a  finite  element  program,  it  is  important 
that  the  convergence  characteri sti cs  of  the  elements  chosen 
for  use  is  known.  To  evaluate  the  convergence 
characteristics  for  the  two  elements,  SH411  and  SH410,  and 
to  determine  the  optimum  mesh  size  and  element  aspect  ratio, 
a  series  of  runs  were  made  on  cylindrical  panels  subject  to 
axial  compression.  The  bifurcation  analysis  with  the  linear 
stress  state  option  was  used  for  all  cases.  The  parameters 
evaluated  were: 

1)  Convergence  character i sties  of  the  SH43  1  and  SH410 
elements. 

2)  Influence  of  panel  ply  orientations  on  the 
convergence  character i sti cs.  Panel  ply  orientations  of 
C45.,-45.32g  and  CO.  ,  45.  ,  -45.  ,  90.  were  evaluated. 

3)  Influence  of  panel  radius  and  panel  size  on  the 
convergence  characteristics. 

4)  Influence  of  mesh  density  and  element  aspect  ratio 
on  the  convergence  characteristics. 

In  Ref  19,  whrn  the  convergence  characteristics  of  the 


two  elements  were  evaluated  for  a  cylindrical  shell,  the 
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SH410  element  converged  very  rapidly  compared  with  the  SH411 
element.  It  was  believed  that  ths  rapid  convergence  was  the 
result  of  case-dependent  cancellation  of  errors  in  different 
directions.  The  SH410  element  also  tended  to  converge  more 
rapidly  than  the  SH411  element  in  this  convergence  study. 
The  reason  for  the  better  performance  of  the  SH410  element 
is  not  evident  although  it  is  probably  related  to  the 
cylindrical  panel  being  curved  in  only  one  direction.  Figure 
11  shows  typical  convergence  results  for  a  C4G.,-45.32S 

t 

panel  using  the  SH411  element  and  element  aspect  ratios  of 
0.5,  1.0,  and  2.0.  Figure  12  compares  the  convergence 
characteristics  of  the  SK411  and  SH410  elements  for  the 
C45. ,-45.32S  and  CO. , 45. , -45. , 90. 3S  panels. 

The  panel’s  ply  orientation  had  the  greatest  influence 
on  the  rate  of  convergence  with  the  C45.,-45.32S  panel 
converging  slower  than  the  £0. , 45. , -45. , 90. 3g  panel.  The 
element  aspect  ratio  also  influenced  the  rate  of 
convergence.  An  element  aspect  ratio  (hei ght/wi dth )  of  1  had 
the  best  rate  of  convergence  followed  closely  by  an  aspect 
ratio  of  0.5.  The  panel  radius  and  panel  size  (with  panel 
aspect  ratio  of  1)  did  not  significantly  influence  the 
convergence  character! sti cs.  Both  the  SH410  and  5H411 
elements  appeared  to  be  converging  monoconical ly  once  an 
element  mesh  size  of  32  elements  was  reached. 

The  accuracy  of  the  STA3S-f31  program  was  evaluated  by 
comparing  the  bifurcation  lead  cf  a  square  CO., 90.  3, g  flat 
plate  with  simple  supported  edges  with  the  theoretical 
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solution  shown  in  Ref  21. 


The  theoretical  critical 

a  e! 

6* 


load  is 


i 


These  results  are  shown  in  Table  IV. 


Table  IV 

Comparison  of  SH410  and  SH411  Elements 
Bifurcation  Loads  to  Theoretical  Results 


Element 

Mesh 

SH410 

Percert 

Error 

SH4 1 1 

Percent 

Error 

6X6 

1.547 

4.74 

1.529 

5.85 

10  X  10 

1.597 

1.79 

1.537 

2.28 

14  X  14 

1.611 

0.80 

1.605 

1.17 

18  X  18 

1.616 

0.49 

1,613 

0.63 

The  SH410  element  required  only  approximately  one-half 
of  the  computer  time  required  for  the  SH411  element.  For 
both  flat  plates  and  the  cylindrical  shells,  the  SH410 
element  has  better  convergence  characteristics  and  a  smaller 
error  than  the  SH411  element. 


For  the  evaluation  of  the  moisture  and  temperature 
influences  on  the  bifurcation  load,  a  mesh  density  of  14 
elements  by  14  elements  was  chosen.  Comparing  the 
bifurcation  load  at  the  14  by  14  element  mesh  with  an 
extrapolated  bifurcation  load  shows  that  for  the 

CO. , 45. , -45. , 90. 1-  panel  the  SH410  element  had  an  error  of 

S3 

0.64  percent  and  the  SH411  element  had  an  error  of  3.38 
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percent.  For  the  E45.,-45.l2g  panel ,  the  SH410  element  had 

an  error  of  10.45  percent  and  the  SH411  element  had  an  error 

of  27.44  percent.  The  extrapolated  bifurcation  load  was 

calculated  using  the  last  three  values  for  the  bifurcation 

a  2 

load  for  each  case  and  then  applying  the  Aitken  A 
"extrapolation  to  the  limit"  method. 

The  SH410  element  will  be  used  in  the  moisture  and 
temperature  evaluation  because  of  its  better  accuracy  and 
quicker  computer  run  times.  The  finite  element  mesh  chosen 
will  be  14  by  14  elements  (196  elements  total)  with  225 
nodes.  The  mesh  is  square  with  an  element  aspect  ratio  of 
1. 


\ 

89 


■’4*. 


Bifurcation  Load,  N  (lbs.) 


Figure  11.  Convergence  Characteristics  of  the  3H411  Element 
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IV..  Moisture  and  Temperature  Conditions  Evaluated 
Moisture  Conditions 

Test  data  (13)  for  the  AS/3501-5  graphite/epoxy  system, 
which  included  the  degradation  in  the  transverse  and  shear 
moduli,  was  available  for  a  saturation  moisture 
concentration  of  1.05  percent.  The  moisture  concentration 
is  measured  as  a  percentage  of  weight  gained.  To  determine 
the  moisture  concentration  through  the  laminate's  thickness 
as  a  function  of  time,  the  series  solution  to  the  Fick 
equation  is  used.  The  coefficients,  and  C 2, 

representing  the  moisture  concentration  at  the  surface  of 

the  laminate,  are  selected  along  with  the  desired 

*  t  2. 

dimensionless  time  T  (where  T  =  K(in  /sec)  x 

2  2 

t(src)/h  (in)  ).  In  the  evaluation  of  the  moisture  and 
temperature  effects,  three  moisture  concentration  conditions 
were  evaluated  at  each  of  five  time  intervals.  These 
moisture  conditions  are  shown  in  Table  V. 

Table  V 

Moisture  Conditions 


Cond.  No. 

co 

C1 

C2 

1 

0.00 

0.00 

0.0105 

2 

0.00 

0.0105 

0.00 

3 

0.00 

0.0105 

0.0105 

The  coefficient  C^  represents  the  initial  moisture 


concentration  in  the  laminate.  For  the  cylindrical  panel , 
Cj  is  the  moisture  concentration  on  the  inside  (-Z)  surface 
and  is  the  moisture  concentration  on  the  outside  (+Z) 
surface.  Conditions  1  and  2  will  result  in  an  unsymmetri c 
degradation  of  the  E^  and  6^  moduli  resulting  in  an 
unsymmetric  laminate.  This  will  introduce  bending-extension 
coupling.  Condition  3  is  symmetric  and  will  not  produce  any 
bending-extension  coupling. 

The  dimensionless  times  used  in  this  analysis  are  0.0, 
0.001,  0.01,  0.1,  and  0.5.  Table  VI  shows  the 
correspondence  between  real  time  and  T*.  The  moisture 
distribution  through  the  thickness  for  moisture  conditions 
1  and  3  and  the  five  time  values  are  shown  in  Fig  13.  A  T* 
of  0.5  represents  the  steady-state  distribution.  The 
moisture  distribution  is  shown  in  Fig  13  as  a  continuous 
function.  As  noted  in  Section  II,  the  moisture 
concentration  is  calculated  at  the  center  of  each  ply  and 
then  assumed  to  be  constant  through  the  thickness. 


Table  VI 


Relation  Between  Real  and  Dimensionless  Time 

Real  Tina  Real  Time  Dimensionless  Time  T* 

(sec)  (days) 


0.0  0.0  0.0 

3.045E04  0.35  0.001 

3.045E05  3.52  0.01 

3.045E06  35.24  0.1 

1.527E07  176.24  0.5 


Note:  Those  times  wore  calculated  using  K  =  0.52537E-10 

2 

(in  /sec)  for  an  3-ply,  0.40  thick,  AS/3501-5  laminate.  The 
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Weight  Gain  (Percent)  Weight  Gain  (Percent) 


Distance  Through  Thickness  (in.) 


Distance  Through  Thickness  (in.) 


Figure  13.  Moisture  Concentration  Distribution  for  Moisture  Conditions 
1  and  3 
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parametric  equation  used  to  determine  K  was  taken  from  Ref 
15.  This  equation  is 


K  =  657  exp(-^jr^)  X  (.03937)Z 

2 

where:  K  -  diffusion  Coefficient  in  in  /sec 
T  -  Temperature  in  degrees  Kelvin 

•l?!DB§!I3iyC§  QBOBi^iSQS 

The  test  data  for  the  AS/3501-5  was  taken  at  four 
temperatures:  80  F,  200  F,  250  F,  and  300  F.  These  four 
temperatures  were  used  to  evaluate  the  influence  of 
temperature.  A  thin  laminate  reaches  temper  ature 

equilibrium  very  quickly  when  compared  with  the  time  to 
reach  moisture  equilibrium.  Thus,  the  laminate  was  assumed 
to  be  at  a  constant  temperature  in  this  evaluation. 


Calculation  of  Moduli  Degradations 

For  a  calculated  moisture  concentration  and  at  a  given 
temperature,  the  reduced  and  G12  values  are  calculated 
from  the  data  shown  in  Table  I.  This  requires  an 
interpolation  of  the  test  data  at  the  given  temperature.  A 
linear  interpolation  between  the  data  at  the  known  moisture 
concentrations  was  used  to  calculate  intermediate  values. 
In  this  manner,  the  and  values  for  each  ply  were 

determined.  For  moisture  condition  3,  the  reduced  moduli 
distribution  through  the  thickness  for  E  and  G  _  are  shown 

jl.  1  4. 
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on  Figs  14a  and  14b,  respectively.  These  two  figures 
represent  the  general  trend  in  moduli  degradation  as  a 
function  of  time  and  temperature  and  do  not  show  the 
degradation  of  moduli  through  the  thickness. 


Laminate  Ply  Orientations 

Three  8-ply  laminates  will  be  evaluated.  The  three  ply 
orientations  chosen  are  CO. ,  45. , -45. , 90.  3^,  C90. ,45. ,- 

45.  ,0.  3S,  and  C45.  ,  -45.  The  relationship  of  these  thr  ee 

to  the  STAGS— Cl  coordinate  system  is  shown  in 


orientations 


STAGS— Cl  Cases 


The  combination  of  moisture  conditions,  temperatures, 
and  times  generates  a  matrix  of  60  cases  for  each  laminate 
per  boundary  condition.  These  60  cases  are  broken  into 
three  sets  of  20  cases;  20  cases  for  each  moisture 
condition.  Because  of  the  large  number  of  cases  to  be 
evaluated,  all  60  cases  were  run  for  only  the  boundary 
condition  with  the  fixed  vertical  edges.  Only  the  20  cases 
for  the  first  moisture  condition  were  run  for  the  vertical 
edge's  simple-supported  boundary  condition  to  evaluate  the 
influence  of  this  change.  The  matrix  of  case  numbers  and 
the  corresponding  conditions  is  shown  in  Table  VII. 


Table  VII 

Moisture  and  Temperature  Conditions  Evaluated 


Case 

No. 

Moi sture 
Cond. 

Boundary 

Cond. 

Laminate 

Element 

Type 

1-20 

1 

Fixed 

CO. ,45. ,-45. ,90. IS 

410 

21-40 

1 

Fixed 

C90. ,45. ,-45. ,0. IS 

410 

41-60 

1 

Fixed 

C45. ,-45. 32S 

410 

101-120 

1 

Simple 

CO.  ,45. ,-45. ,90. 3S 

410 

121-140 

1 

Simple 

C90. ,45. ,-45. ,0. 3S 

410 

141-160 

1 

Simple 

C45. ,-45. 32S 

410 

201-220 

2 

Fixed 

CO. ,45. ,-45. ,90. 3S 

410 

221-240 

2 

Fixed 

C90. ,45. ,-45. ,0. 3S 

410 

240-260 

2 

Fixed 

C45. ,-45. 32S 

410 

401-420 

3 

Fixed 

CO. ,45. ,-45. ,90. 3S 

410 

420-440 

3 

Fixed 

C90. ,45. ,-45. ,0. 3S 

410 

441-460 

3 

Fixed 

C45. ,-45. 32C 

410 

50 


These  cases  were  all  run  -for  12  in.  by  12  in.  panels 
with  a  panel  radius  of  12  in.  Two  additional  sets, 
identical  to  cases  1-20,  were  run  with  panel  radii  of  24  in. 
and  48  in.  to  evaluate  the  influence  of  panel  radius  on  the 
degradation  characteri sti cs. 

The  computation  of  the  moisture  concentration 
distribution  and  the  reduction  in  moduli  was  done  by  a 
computer  program.  A  listing  of  this  computer  program  is 
included  in  Appendix  A.  This  computer  program  also 
generated  the  input  deck  for  the  STA6S-C1  program. 

Once  the  value  of  was  calculated,  the  corresponding 
value  of  was  calculated  using  the  relationship 

v,z  _  Vz, 

f(  ~  £2 

As  noted  previously,  both  E1  and  V"12  were  assumed  to  be 
constant.  Actually,  both  of  these  material  properties  vary 
with  moisture  concentrati on  and  temperature.  E,  increases 
by  about  10  percent  and  ^12  varies  from  0.34  to  0.45. 
Using  the  STAGS-C1  program,  the  increase  in  E^  resulted  in 
an  approximately  eight  percent  increase  in  the  bifurcation 
load.  l/*j2  did  not  significantly  influence  the  bifurcation 
load  with  changes  of  less  than  two  percent. 

Since  caused  an  increase  in  the  bifurcation  load, 
it  was  decided  not  to  include  this  change  in  the  case 
studies  because  it  would  make  an  evaluation  of  the  effects 
of  reductions  in  the  other  two  moduli  more  difficult. 
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V._  Results  and  Discussion 

Reduction  in  Bifurcation  Load 

The  moisture-  and  temper ature-i nduced  degradations  in 

the  E2  and  G,  moduli  resulted  in  reductions  in  the  panels’ 

bifurcation  load.  The  results  of  the  STAGS-C1  runs  for  the 

three  panels  with  the  fixed  boundary  condition  are  shown  in 

Figs  16,  17,  and  18.  In  these  plots  Nx  represents  the 

on  g 

bifurcation  load  for  the  room  temperature  condition  at 
T=O.0.  This  condition  is  unaffected  by  either  temperature- 
or  moisture  degradations.  These  results  are  also  shown  in 
tabular  form,  along  with  individual  plots  fo-  each  time  end 
temperature  series,  in  Appendix  B.  These  results  are  for  12 
in.  by  12  in.  panels  with  a  radius  of  12  in. 

As  was  expected,  the  panels’  bifurcation  load  decreased 
with  increasing  temperature  and  absorbed  moisture.  At  the 
higher  temperatures  and  moisture  concentrations,  this 
reduction  was  significant;  ranging  from  21  percent  for  the 
CO. , 45. , -45. , 90. 3g  laminate  to  43  percent  for  the  C45. , 
-45.32S  laminate.  These  reductions  are  especially 
significant  considering  that  the  longitudinal  moduli  was  not 
changed.  A  summary  of  the  maximum  reduction  in  for  each 
laminate  and  moisture  condition  is  shown  in  Table  VIII. 
These  results  are  for  the  fixed  boundary  condition. 


Figure  18.  Degradation  in  Nx  for  the  (^5* )2g  Laminate 
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Table  VIII 


Percent  Reduction  in  Bifurcation  Load 
at  300  F  and  T  =0.5 

Moisture  Condition 

Laminate  1  23 

CO. ,45. ,-45. ,90. IS  13.4  15.5  21.3 

C90. ,45. ,-45. ,0. IS  15.9  17.7  24.5 

C45. , —45. I2S  20.6  29.9  42.7 

The  change  in  boundary  conditions  did  not  significantly 
change  the  panels’  bifurcation  loads.  Comparisons  of  the 
calculated  N^s  for  the  two  boundary  conditions  are  shown  in 
Table  IX.  These  results  are  shown  for  the  original  room 
temperature  values  and  the  worst  case  values  at  300  F  and  T* 
=0.5  for  moisture  condition  1. 

Table  IX 


Change  in  N  Due  to  Changing  Boundary  Conditions 
X 

Moisture  Cond.  1,  T*  =  0.0,  80  F 


Laminate 

Boundary 

Fixed 

Condition 

Simple 

Percent 

Reduction 

CO. ,45. ,-45. ,90. IS 

514.8 

502.8 

-2.33 

C90. ,45. ,-45. ,0. IS 

446.0 

444. 1 

-0.43 

C45. ,-45. 12S 

428.9 

421.6 

-1.63 

Moisture  Cond.  1,  T*  =  0.5, 

300  F 

Lami nate 

Boundary 

Fixed 

Cond4.  ti  on 
Simple 

Percent 

Reduction 

CO. ,45. ,-45. ,90. IS 

445.9 

434.6 

-2.53 

C90. ,45, , -45. ,0. IS 

375.3 

373.5 

-0.48 

C45. ,-45. 12S 

340.4 

328.9 

-3.78 
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Moisture  conditions  1  and  2  cause  the  initially 
symmetric  laminate  to  become  unsymmetric.  This  unsymmetry 
introduces  bending-extension  coupling.  The  influence  of 
this  coupling  can  be  seen  in  Figs  16,  17,  and  18.  For  the 

CO. , 45. , -45. , 90. Dg  laminate  and  the  C90. , 45. , -45. , 0. Dg 
laminate,  this  coupling  does  not  significantly  change  the 
panels’  .  There  is  only  a  three  percent  difference 
between  the  two  moisture  conditions  for  these  two  laminates. 
The  difference  between  the  two  conditions  for  the  C45. , 
-45. D^g  laminate  is  10  percent  at  300  F  ana  a  time  of  0.5. . 
The  influence  of  the  coupling  stiffnesses  is  greater  for 
this  laminate  because  the  magnitude  of  the  terms  are 

larger  than  they  were  for  the  other  two  laminates  while  the 
magnitude  of  this  laminate’s  and  Di  j  terms  are  generally 

smal 1 er . 

To  evaluate  the  influence  of  panel  radius  on  the 
moisture-  and  temperature-caused  reductions,  the  20  cases 
for  moisture  condition  1  for  the  CO. , 45. , -45. , 90. 3g  panel 
were  rerun  with  radii  of  24  in.  (Cases  601-620)  and  48  in. 
(Cases  701-720).  As  expected,  the  panel’s  original 

reduced  with  increasing  radii.  changed  from  514.8  lbs. 

at  a  12  in.  radius  to  290.8  lbs.  at  a  24  in.  radius  to  170.3 
lbs.  at  a  46  in.  radius.  The  percentage  of  reduction  from 
this  original  did  not  significantly  vary  from  those  for 
the  12  in.  radius  panel .  This  indicates  that  the  moisture- 
anri  temperature-caused  reductions  are  not  significantly 
influenced  by  panel  radius;  hence,  the  results  obtained  for 
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the  12  in.  radius  panel  should  bs  valid  -for  any  radius 
within  the  limits  of  a  shallow  shell. 

Pre-Buckl ed  Displacement  Character  i_oti  cs 

Unlike  flat  plates  and  straight  columns,  the 
cylindrical  shell  experiences  normal  di spl acements  and 
rotations  prior  to  buckling.  The  pre-buckled  displacements 
for  the  three  panels  were  very  similar  with  the  differences 
being  one  of  magnitude.  For  the  CO. , 45. , -45. , 90. 1^  and 
C90. , 45. , -45. , 0. panels,  the  di spl acements  were  small  and. 
within  the  limits  of  the  small  displacement  assumptions. 
The  displacements  for  the  C45.  , -45.  panel  were  an  order 

of  magnitude  larger.  A  typical  plot  of  the  pre-buckled 

displacement  geometry  for  the  CO. , 45. , -45. , 90. panel  is 
shown  in  Fig  19.  This  is  a  plot  for  Case  No.  20  of  the  U, 
V,  W  resultant  component  displacements.  The  scale  factor  is 
20,000.  Contour  plots  representing  lines  of  constant 
displacement  or  rotations  for  the  U,  V,  W,  RU,  and  RV 

displacements  and  rotations  for  this  case  are  shown  in  Figs 
20-24.  The  numbers  adjacent  to  the  contour  lines  signify 
the  percentage  of  maximum  displacement  that  that  contour 
line  represents.  The  influence  of  the  fixed  boundary 
conditions  are  clearly  shown  in  the  contour  plots. 

Additional  displacement  and  contour  plots  for  representati ve 
cases  for  each  ply  orientation  and  boundary  condition  are 
included  in  Appendix  C. 

The  ply  orientation  and  boundary  conditions  controlled 
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the  overall  shape  of  the  pre-buckled  displacement.  The 
moduli  degradations  due  to  moisture  and  temperature  did  not 
change  the  displacement  shape,  but  did  influence  the 
magnitude  of  the  displacements.  Figure  25  shows  the  W 
displacements  along  row  8,  which  is  in  the  circumf erential 
direction.  Comparison  plots  for  several  cases  representing 
both  fixed  and  simple  boundary  conditions  for  the  CO. ,45., 
-45., 90. laminate  are  shown.  Similar  plots  for  the 
C90. , 45. , -45. , 0. laminate  are  shown  in  Fig  26. 
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CRSE  NO.  20 

PLOT  NO.  1.  UNIT  0.  STEP  0 
DISPLACEMENT  GEOMETRY 
UVW  COMPONENT 

MODEL  SCRLE  =  .5000E*00.  ORIENT.  =  0.00.  60.00. 

SOLUTION  SCRLE  =  .20Q0E*05 


60 


CASE  NO.  20 

PLOT  NO.  3,  UNIT  1,  STEP  0 
DISPLACEMENT  CONTOURS 
U  COMPONENT 

HODEL  SCALE  =  .5000E+00.  ORIENT.  =  0.00.  0.00. 

SOLUTION  SCALE  =  . 2273E+07 
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CASE  NO.  20 

PLOT  NO.  5 ,  UNIT  1.  STEP 
DISPLACEMENT  CONTOURS 
M  COMPONENT 

MODEL  SCALE  =  .SOOOE+OO,  ORIENT.  = 
SOLUTION  SCALE  =  .4701E+07 


Figure  22.  W  Component  Mspi?  cement  Contour  Plot,  for  (0.  »-^5»  *90.  )£ 
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CRSE  NO.  20 

PLOT  NO.  6.  UNIT  1 .  STEP  0 
DISPLACEMENT  CONTOURS 
RU  COMPONENT 

MODEL  SCALE  =  .5000E+00,  ORIENT.  =  0.00.  0.00 

SOLUTION  SCALE  =  . 47S7E+Q7 
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CRSE  NO.  20 
PLOT  NO.  7.  UN 
DISPLRCEMENT  CO 
RV  COMPONENT 

HOOEL  SCALE  =  .SOOOE 
SOLUTION  SCALE  =  .47 


Fixed  Boundary  Condition 


Cato  No.  1 

T*  =  0.0 
v  80  F 
^M.C.  1 


Case  No.  20 
* 

T  ■=  0.5 
300  F 

M.C.  1 


Case  No.  d20 


T  =  0.5 
300  F 
M.C.  3 


Simple  Boundary  Condition 


Case  No.  101 


T  =  0.0 
80  F 

M.C.  1 


Notei  Displacements  are  of  magnitude  10  (in.)  for  a  unit 

axial  load;  M.C.  is  the  moisture  condition. 


Figure  25.  Comparison  of  Pre-buckled  W  Displacements  for  (0. ,45. ,-45. ,90. ) 


Laminate 


(0. ,45. ,-45. ,50. )_  Fixed  Boundary  Condition 


Case  No.  ] 

* 

T  =0.0 
cO  F 

V  p  1 
i'«  ,  L.  • 


(90. ,45. ,-45. »0.)_  Fixed  Boundary  Condition 


Case  No.  21 

T*  »  0.0 
80  F 
K.C.  1 


(90.  ,45.  ,-45.  ,0.  )c.  Simple  Boundary  Conditioi 


Case  No.  121 
T*  =  0.0 

so  f 

*  M.C.  1 


(90. ,43. , -45. ,0. )„  Fixed  Boundary  Condition 


Case  No.  440 


T  =  0.5 
300  F 

M.C.  3 


-4  t  \ 

Note:  Displacements  are  of  magnitude  1C  (in.;  for  a  unit  axial 
load;  M.C.  is  the  moisture  condition. 


Figure  26.  Comparison  cf  Fre-buckled  W  Displacements  for 

{0. ,45.. -45., 90.),  and  (90. ,45, ,-45* »0. )s  Laminates 
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In  comparing  the  plots  for  the  two  laminates,  it  is 
seen  that  the  CO. , 45. , -45. , 90. 3g  laminate  develops  a 
moderate  but  significant  hump  along  the  vertical  edge  for 
the  fixed  boundary  condition.  In  Table  X,  the  .  bending 
stiffnesses  for  these  two  laminates  are  shown. 

Table  X 

Bending  Stiffnesses  for  the  CO. , 45. , -45. , 90. 3^ 
and  C90. , 45. , -45. , 0. 3g  Laminates  at  80  F  and  T*  »  0.0 


CO. ,45. ,-43. ,90. 3g 

C90. , 45. , -45. , ' 

11 

0. 72001E02 

0. 19337E02 

12 

0. 11012E02 

0. 1 1012E02 

16 

0. 43S87E01 

0. 43B87E01 

22 

0. 19337E02 

0. 7200) E02 

26 

0.43837E01 

0. 43887E01 

66 

0. 13297E02 

0. 13297E02 

The  only  difference  between  the  two  laminates  is  the 

Djj  and  D^.,  values.  Since  the  ^  values  are  the  same  and 

the  B.  .  values  are  zero  for  this  case,  it  appears  that  the 
i  j 

larger  hump  in  the  CO. , 45. , -45. , 90. 3g  laminate  is  because 
its  ®22  stiffness  is  smaller  than  the  C90. , 45. , -45. , 0. 3g 
laminate.  Tt  is  possible  ths>t  these  two  humps  may  be 

acting  as  two  connected  cylindrical  panels  with  a  smaller 
•  ;  radius  which,  along  with  the  larger  stiffness 

.  !  iv.nstf',  provides  the  additional 


'  L  -• 


the 


CO. . 45. , -A5. , 90. 3 


laminate  to  have  a  higher  bifurcation  load. 

The  simple  supported  boundary  condition  does  not 
provide  rotational  restraint  along  the  vertical  edges.  The 
plots  for  Case  No.  101  and  121  in  Figs  25  and  26  show  these 
boundary  conditions.  The  simple  supported  edge  reduces  the 
size  of  the  hump  in  each  case  but  the  ICO.  ,  45.  ,  -45.  ,  VO. 
laminate  appears  to  have  a  slightly  smaller  effective  radii 
which  results  in  a  higher  bifurcation  load. 

The  shape  of  the  pre-buckled  displacement  for  the 
C45.,-45.12g  laminate  is  similar  to  those  for  the  other  two 
laminates.  However,  the  displacements  are  an  order  of 
magnitude  larger.  Figure  27  shows  the  displacement  plot  of 
the  U,  V,  W  displacement  component  at  the  same  scale  factor 
as  the  previous  plot.  Figure  28  shows  the  plots  of  the  W 
displacement  in  the  circumf erenti al  direction.  The  drastic 
change  in  magnitude  of  the  W  displacement  is  clearly  seen  by 
comparing  Case  No.  460  with  Case  No.  41. 

In  order  to  determine  which  of  the  reductions  in  moduli 
is  causing  the  large  deflections  for  the  £45. , -45. 

a.u 

laminate,  two  STAGS-C1  runs  were  made  in  which,  separately, 
the  E2  and  moduli  were  reduced  to  20  percent  of  the 

original  values.  Plots  for  the  W  displacement  in  the 
circumferential  direction  for  both  fixed  and  simple  boundary 
conditions  are  shown  in  Fig  29.  The  reduction  in  the 
moduli  causes  significantly  larger  displacements  than  a 
reduction  in  F_,  does.  This  i-auid  increase  in  displacement 
with  reduced  moduli  explains  the  large  reductions  in 


Fixed  Boundary  Condition 


_o 

Notes  Displacements  are  of  magnitude  10  J  (in.)  for  a  unit 
axial  load. 


Figure  29. 


Comparison  of  Pre-Buckled  W  Displacements  with  E^  and 
Reduced  to  20  Percent  for  the  (45,  ,-45.  )^g  Laminate 


noted  above 


The  maximum  displacements  -for  the  C45.,-45. 3^g 

laminate,  at  the  bifurcation  load,  are  of  the  order  of  10  1 

-2 

with  extensional  strains  of  the  order  of  10  and  rotations 
approaching  5  degrees.  These  large  displacements  are 
approaching  the  limits  of  small  displacement  theory 
associated  with  the  choice  of  the  bifurcation  analysis  used 
in  this  analysis  displacement.  A  more  rigorous  analysis 
using  the  STAGS-C  finite  difference  program  with  a  non¬ 
linear  buckling  analysis  would  be  appropriate  for  this 
laminate. 

Eigenvector 

Each  of  the  three  laminates  has  a  distinctive 
eigenvector.  Displacement  and  contour  plots  for  each 
laminate’s  eigenvector  are  shown  in  Figs  30-35.  The 
amplitude  of  the  eigenvector  in  the  displacement  plots  has 
been  exagerated  to  shew  the  shape. 

The  CO. , 45. , -45. , 90. lg  laminate  develops  a  series  of 
waves  that  are  oriented  with  the  direction  of  the  outer  45 
degree  plies.  The  amplitude  of  the  waves  increases  to  a 
maximum  near  the  center  of  the  panel.  Moisture  and 
temperature  changes  and  the  boundary  conditions  do  not  alter 
the  shape  of  the  eigenvector.  The  C90. , 45. , -45. , 0. 3g  and 
C45.  ,-45.Il2g  laminates  develop  a  series  of  axial  waves. 
Again,  neither  the  temperature  and  mos store  changes  or  the 
boundary  conditions  alter  the  shape  of  the  eigenvector. 


CASE  NO.  20 

PLOT  NO.  2.  UNIT  0,  STEP  0 
EIGENVECTOR  GEOMETRY 
UVW  COMPONENT  .  MODE  1 

MODEL  SCALE  =  .5000E+0Q.  ORIENT.  =  0.00,  so  00 

SOLUTION  SCALE  =  .5000E+00 


Figure  30.  Eigenvector  Geometry  for  (0. ,45. ,-45. ,90. )  Laminate 


CASE  NO.  20 

PLOT  NO.  10.  UNIT  1.  STEP  0 
EIGENVECTOR  CONTOURS 
N  COMPONENT  ,  MODE  1 

MODEL  SCALE  =  .SOOOE+OO.  ORIENT.  =  0.00.  0.00. 

SOLUTION  SCALE  =  .3497E+02 


CASE  NO.  21 

PLOT  NO.  2.  UNIT  0.  STEP  0 
EIGENVECTOR  GEOMETRY 
UVN  COMPONENT  ,  MODE  1 

M0OEL  SCALE  =  .5000E+00,  ORIENT.  =  0.00.  60.00. 

SOLUTION  SCALE  =  .5000E+00 
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CASE  NO.  21 

PLOT  NO.  10,  UNIT  1.  STEP  0 
EIGENVECTOR  CONTOURS 
W  COMPONENT  ,  MODE  1 

HOOEL  SCALE  =  .5000E»00.  ORIENT.  =  0.00.  0.00. 

SOLUTION  SCALE  =  . 1529E,03 
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CRSE  NO .  41 

PLOT  NO.  2,  UNIT  0,  STEP  0 

EIGENVECTOR  GEOMETRY 

UVU  COMPONENT  ,  MODE  1 

riQDEL  SCALE  -  .  bOOCE  +  CC.  CLIENT.  G.G 

SOLUTIGN  SCALE  -  . bCQGErCG 


Figure  34.  Eigenvector  Geometry  for  (45.  »-45*  Laminat 


CASE  NO .  . 

PLOT  NO.  10.  UNIT  1 .  STEP  0 
EIGENVECTOR  CONTOURS 
U  COMPONENT  ,  MODE  1 

FCGEX  cC^LE  -  .  SCuO£*CG.  CLIENT.  --  O.CC 

sol u r : cm  sr.qit  -  . ss he +02 


Linear  Approximation  of  Reductions  i_n  the  Bifurcation  Load 

The  W  and  U  displacements  fo''-  the  CO.  ,  45.  ,  -45.  ,  90. 
and  C90. , 45. , -45. , 0. Ig  laminates  were  found  to  vary 
approximately  linearly  with  the  changes  in  due  to 

moisture  and  temperature.  Figure  36  shows  plots  of  W, 

measured  at  the  middle  of  the  panel,  and  U,  measured  at  the 
middle  of  the  top  edge,  vs  the  N#  for  the  CO. , 45. , -45. , 90. 

X  o 

laminate.  The  closed  symbols  represent  the  fixed  boundary 
conditions  and  the  open  symbols  represent  the  simple 
boundary  conditions.  The  lines  connect  the  symmetric  cases, 
represented  by  circles,  while  the  triangle  represents  the 
unsymmetric  moisture  condition  1  and  the  square  represents 
moisture  condition  2. 

The  linear  behavior  of  these  plots  indicated  that  the 
changes  in  due  to  moisture  and  temperature  could  be 

approximated  by  a  linear  function.  It  was  noted  previously 
that  several  additional  STAGS-C1  runs  were  made  with  the  E^ 
and  G,_  moduli  individually  reduced  for  the  C45.  ,-45.3„ 
laminate.  Similar  runs  were  also  made  for  the  C0.,45. ,- 
45.,90.1S  and  C90. , 45. , - 45. , 0. 3g  laminates.  A  plot  of  the 
results,  Fig  37,  showed  that  was  almost  a  linear  function 
of  the  reductions  in  moduli.  Linear  relationships  between 
the  reduced  and  moduli  and  were  developed.  In 

these  equations,  the  is  calculated  as  a  function  of  the 
reduced  moduli  which  is  represented  as  a  percentage  of  the 
original  moduli  value.  3y  combining  these  two  linear 
functions,  one  equation  relating  the  to  the  reduced 
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Bifurcation  Load.  N  (lbs.)  Bifurcation 


A  Moisture  Condition  1,  Fixed 

B  Moisture  Condition  2,  Fixed 

•  Moisture  Condition  3>  Fixed 

A  Moisture  Condition  1,  Sir, pie 

O  Moisture  Condition  3>  Sir.ple 


Figure  37.  Influence  of  E0  and  G10  Degradation  on  N  for 

c  . Yc.  ~  x 

(0.  ,45. ,  -43.  »90.  )c  La-ninate 


moduli  was  developed.  This  equation  is  dependent  upon  the 
ply  orientations  and  boundary  conditions,  but  independent  of 
the  time,  temperature,  moisture  concentrati ons,  or  panel 
radius  within  the  limits  discussed  above. 

To  calculate  the  at  a  particular  time  and 

temperature,  the  Fick  equation  is  used  to  calculate  •  the 
moisture  concentration  distribution  through  the  thickness. 
Using  these  values  and  the  temperature,  the  reduced  moduli 
are  determined  -from  the  test  data.  An  average  reduced 
moduli  value,  represented  as  a  percentage  o-f  the  original 
moduli,  is  calculated  -for  the  laminate  for  both  and  G^* 
These  values  are  input  into  the  equation  for  the  particular 
laminate  ply  orientations  and  boundary  conditions.  The 
panel's  N^,  as  a  percentage  of  the  original  N^,  is 
calculated. 

The  three  ply  orientations  and  two  boundary  conditions 
used  in  this  thesis  require  six  equations  to  predict  the  . 
The  general  form  of  the  equation  is  shown  below  along  with  a 
table  of  the  slope  and  intercept  values  for  each  of  the  six 
cases. 

\  =  ml  E2  +  bl  +  m2  G12  +  b2  ”  100 

where*  N  -  reduced  bifurcation  load  expressed  as 
percentage  of  the  original  R^ 

E„  -  reduced  transverse  moduli  expressed  as 
percentage  of  the  original 

G,,,  -  reduced  shear  moduli  expressed  as  a 

percentage  of  the  original 


For  the  reasons  discussed  previously,  this  equation 


does  not  include  the  effects  of  any  change  in  E^  or 

However,  the  equation  could  be 

modified  to  include 

these 

changes. 

Tabl  e 

XI 

Slope  and  Intercept  Values 

Laminate  E2 

G12 

ml 

bl 

m2 

b2 

CO. ,45. ,-45. ,90. 3S 

Fixed  0.1532 

84. 68 

0. 1094 

89.06 

Simple  0. 1554 

84.46 

0. 1099 

89.01 

C90. ,45, ,-45. ,0. Dg 

Fixed  0.2099 

79.01 

0. 1031 

E9.69 

Simple  0.2097 

79.03 

0. 1044 

89. 56 

C45. ,-45. 12S 

Fixed  0.0772 

92.28 

0.4616 

53.84 

Simple  0.0759 

92.41 

0.4776 

52.24 

The  values  of  and  m0  allow  comparisons  to  be  made  of 
how  equal  changes  in  the  values  of  E 2  and  GJ2  influence  the 
change  in  N^.  It  was  noted  above  that  the  C90.,45., 

-45. ,0. laminate  loses  more  buckling  strength  than  the 
CO. , 45. , -45. , 90. 3g  laminate  does  for  each  moisture  and 
temperature  conditions.  From  the  above  values  of  and  m2 
for  these  two  laminates,  the  apoarent  reason  for  this  is 
that  the  m.  for  the  C90, , 45. , -45. , 0. 3  laminate  is  about  30 

<1  3 

percent  larger  than  the  m.  value  for  the  CO. , 45. , -45. , 90. 1_ 

■4.  3 


f 


laminate.  Thus  -for  equal  changes  in  E^,  the  C90.,45., 
-45.  ,0.1g  laminate’s  bi-furcation  load  degrades  more  rapidly. 

The  comparison  between  the  m^  and  m^  values  for  the 
C45. , —45. 3oc  laminate  supports  the  conclusion  reached  above 
that  the  change  in  GJ2  is  the  major  cause  of  the  large 
reductions  in  this  panel’s  . 

Comparisons  of  the  values  calculated  using  this 
equation  and  the  STAGS-C1  calculated  Ny  values  for  Cases  1- 
20  are  shown  below.  Additional  tables  of  comparisons  for 
other  cases  are  shown  in  Appendix  D.  In  most  cases,  the 
maximum  error  was  less  than  5  percent.  However,  in  some  of 
the  C45. , -45. cases,  the  error  reached  10  percent.  This 
may  be  because,  at  the  higher  temperatures  and  moisture 
conditions,  the  resulting  large  displacements  and  rotations 
are  influencing  the  accuracy  of  the  STAGS-C1  program. 
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Table  XII 

Comparison  of  Predicted  and  Calculated  Values 


Case 

STAGS-C1 

Predicted 

Percent 

No. 

N  Value 

X 

N  Value 

X 

Error 

1 

514.8 

514.80 

.00 

2 

514.2 

514.54 

.07 

3 

512.5 

513.59 

.21 

4 

510.5 

510.99 

.  10 

5 

509.7 

509.57 

-.03 

6 

493.3 

492.23 

-.22 

7 

492.  1 

491.64 

-.09 

8 

488.7 

489.34 

.  13 

9 

484.5 

483.07 

-.30 

10 

482.3 

479.47 

-.59 

11 

482.3 

481.10 

-.25 

12 

479.8 

479.78 

.00 

13 

472.3 

474.71 

.51 

14 

462.8 

460.86 

-.42 

15 

458.0 

452.93 

-1.11 

16 

480,4 

479.19 

-.25 

17 

477.  1 

477.35 

.05 

18 

466.8 

470.05 

.70 

19 

453.2 

450.20 

—  .  66 

20 

445.9 

438.59 

-1.64 
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VK  Conclusions 

On  the  basis  o-f  the  finite  element  analysis,  the 
following  conclusions  can  be  made  for  cylindrical,  composite 
panels  subject  to  moisture  exposure  and  elevated 
temperatures  and  loaded  in  axial  compression. 

1.  The  results  of  the  STAGS-C1  finite  element 
analysis,  using  modified  flat  plate  elements,  compare  well 
with  the  STAGS-C  finite  difference  analysis,  using  a  curved 
surface  representation,  in  predicting  the  bifurcation  load 
of  cylindrical  panels  using  the  linear  pre-buckled 
displacement  analysis  mode. 

2.  The  simpler  4  noded  quadrilateral  finite  element 
<SH410)  had  better  convergence  characteri sties,  was  more 
economical  to  use,  and  yielded  results  comparable  to  the  4 
noded  quadri 1 ateral  with  4  mid-side  nodes  (SH411)  when 
calculating  the  panel's  bifurcation  load. 

3.  The  bifurcation  load  of  a  composite  panel,  with  a 
resin  material  whose  elastic  moduli  are  reduced  by  absorbed 
moisture  and  elevated  temperatures,  will  degrade  with 
increasing  moisture  concentrations  and  temperatures. 

4.  The  extent  of  the  degradation  in  the  bifurcation 
load  is  influenced  by  the  degree  of  moisture  concentration, 
the  temperature,  and  the  panel's  ply  orientations.  At  300  F 
and  a  moisture  weight  gain  of  1.05  percent,  the  CO. ,45. , 
-45. ,90.!g  panel  experienced  a  21.3  percent  degradation,  the 
C90. , 45. , -45. , 0. D_  panel  experienced  a  24,5  percent 

b 

degradation,  and  the  C45. , -43 . panel  experienced  a  42.7 
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w 

\ 

percent  degradation. 

5.  A  change  in  the  rotational  restraints  of  the 
vertical,  straight  sides,  from  fully  fixed  to  simple- 
supported,  did  not  significantly  reduce  the  bifurcation  load 
or  change  the  moisture-  and  temperature-induced  degradation 
characteristics  for  a  panel  aspect  ratio  of  1 . 

6.  Increasing  the  cylindrical  panel's  radius  decreased 
the  panel’s  bifurcation  load  but  did  not  significantly 
change  the  moisture-  and  temperature-induced  degradation 
characteristics. 

7.  The  bending-extension  coupling  resulting  from  the 
unsymmetric  moisture  concentration  distributions,  which 
resulted  in  an  unsymmetric  laminate,  did  not  significantly 
influence  the  bifurcation  load. 

8.  The  cylindrical  panel’s  bifurcation  load  is 
influenced  by  the  ply  orientations  and  the  panel’s  pre¬ 
buckled  displacements  which  vary  with  different  ply 
orientations. 

9.  The  reduction  in  the  bifurcation  load  due  to 

moisture  and  temperature  was,  in  general,  linearly  related 
to  the  average  reduction  in  the  transverse,  and  shear, 

moduli.  A  linear  relationship  can  be  developed  using 
the  results  of  three  STAGS3-C1  runs  for  each  ply  orientation, 
panel  aspect  ratio,  and  boundary  conditions  from  which  the 
bifurcation  load  can  be  calculated  from  the  average  E,,  and 
*  Gj^  values  for  any  moisture  and  temperature  condition. 
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Agger tlix  B 

Summary  of  STASE-Cf  Runs 


Case  No.  1-20  CO. , 45. , -45. , 90. 3S 

Laminate  Nondimensional  Time 

Temperature 


<Deg.  F> 

!  0.00  ! 

0.001  5 

o.oi  : 

0.  1  ! 

0.5 

80.0 

514.8 

514.2 

512.5 

510.5 

509.7 

1 

.9988345 

.9955322 

.9916472 

. 9900932 

200.0 

493.3 

492.  1 

48B.7 

484.5 

482.3 

. 9582362 

. 9559052 

.9493007 

.9411422 

.9368687 

250.0 

482.3 

479.8 

472.3 

462.8 

458 

. 9368687 

.9320124 

.9174437 

. 8989399 

. 8396659 

300.0 

480.4 

477.  1 

466.8 

453.2 

445.9 

.9331779 

.9267677 

.9067599 

.8803419 

.8661616 

Case  No. 

Laminate 
Temperature 
(Deg.  F> !  0.00 

21-40  C90. , 45. , — 

Nondimensional  T 

S  0.001  !  0.01 

45. ,0. 3S 

ime 

:  0. 1  ! 

0.5 

80.0 

446 

445.  1 

442.9 

440.9 

440.  1 

1 

.9979321 

. 9930493 

.9885650 

.9867713 

200.0 

422.2 

420.9 

417 

412.9 

410.6 

. 9466368 

.9434978 

.9349776 

.9257348 

. 7206278 

250.0 

412 

409 

400.5 

392.  1 

386.  7 

. 9237668 

.9170404 

.8979821 

.8791480 

. 8670404 

300.0 

410.5 

406.3 

395.4 

382.7 

375.3 

.9204036 

.9121076 

.8865471 

.3580717 

.8414798 

Laminate 
Temperature 
(Deg.  F)  i 

Case  No.  41-60  C45.,-4532S 

Nondimensional  Time 

o.oo  :  o.ooi  i  o.oi  :  o.i  : 

0.5 

80.0 

428.9 

428.6 

427.7 

426.7 

426.4 

1 

. 9993005 

.9972021 

.9948706 

.9941711 

200.0 

408.  4 

407 . 8 

406.  1 

401.4 

397.7 

• 

9522033 

.9508044 

.9463*08 

. 9358825 

. 9272558 

250.0 

337.4 

336.  3 

332.  5 

371.4 

362.2 

a 

9032403 

, 9006761 

.8918163 

.9659361 

. 3444939 

300.0 

381.7 

380 . 3 

375 

356.6 

340.  4 

■ 

8899510 

. 8666869 

. 8743297 

.8314292 

. 7936532 
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Case  No.  101-120  CO. , 45. , -45. , 90. 3S 


Laminate  Nondimensional  Time 

Temperature 


(Oeg.  F) 

:  o.oo  i 

0.001  : 

0.01 

o.  i  : 

0.5 

eo.o 

502.8 

502.2 

500.5 

498.5 

497.6 

1 

.9988067 

.9954256 

.9914479 

. 9896579 

200.0 

481.2 

480 

476.6 

472.4 

470.3 

. 9570406 

.9546539 

.9478918 

.9395336 

. 9353620 

250.0 

470.3 

467.8 

460.3 

451 

446.3 

.9353620 

. 9303898 

.9154733 

.8969769 

. 8876293 

300.0 

468.5 

465.2 

455 

441.7 

434.6 

.9317820 

.9252188 

.9049324 

. 8784805 

. 8643596 

Case  No.  121-140  C90. , 45. , -45. , 0. ]S 

Laminate 

Nondimen 

sicnal  Time 

Temperature 

(Deg.  F) 

!  0.00  5  C.001  ! 

0.01  !  0.1  ! 

0.5 

80.0 

444.1  443.3 

441.2  439.2 

438.4 

1  .9981986 

.9934699  .9889664 

.9871651 

200.0 

420.6  419.3 

415.5  411.5 

409.2 

.9470840  .9441567 

.9356001  .9265931 

.9214141 

250.0 

410.5  407.7 

399  339.9 

385 

.9243414  .9180365 

.8934463  .0779554 

.8669219 

300.0 

409  405.5 

393.8  380.9 

373.5 

.9209637  .9130326 

.8867372  .8576897 

.8410268 

Case  No 

Laminate 
Temperature 
(Deg.  F) !  0.00  ! 

.  141-160  C45. , —45. 32S 

Nondimensional  Time 

0.001  !  0.01  !  0.1 

0.5 

80.0 

421.6 

421.3 

420.5 

419.6 

419.2 

1 

.9992884  . 

9973909  . 

9952562 

.9943074 

200.0 

401 

400.5 

398.6 

393.6 

389.8 

.9511385 

.9499526  . 

9454459  . 

9335363 

.9245731 

250.0 

379.7 

378.5 

374.3 

362.3 

352.9 

.9006167 

.8977704  . 

8373093  . 

3593454 

. 8370493 

300.0 

373.9 

372.3 

366 

345.  8 

328 . 9 

.  8868576 

.8330645  . 

8631214  . 

3202007 

. 7801233 
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Case  No.  201-220  [0. , 45. , -45. , 90. IS 

Laminate  Nondi mensi onal  Time 

Temperature 


(Deg.  F) 

!  0.00  5  0.001  ! 

0.01  : 

0.  1 

0.5 

80.0 

514. B  514.4 

513.3 

511.7 

510.7 

1  .9992230 

.9970862 

.9939782 

. 9920357 

200.0 

493.3  491.9 

487.5 

482.7 

480.  3 

.9582362  .9555167 

.9469697 

.9376457 

.9329837 

250.0 

482.3  479.2 

468.9 

458 

452.6 

.9360687  .9308469 

.91 0S392 

.8896659 

.8791764 

300.0 

480.4  475.6 

460 

443.  1 

435.  1 

.9331779  .9238539 

. 8935509 

.8607226 

.8451826 

Case  No.  221-240  C90. , 45. , -45. , 0. US 

Laminate  Nondi mensi onal  Time 

Temperature 

(Deg.  F) i  0.00  1  0.001  !  0.01  !0.1  ! 

in 

« 

o 

80. 0 

445.9 

445.  1 

443.2  441.9 

440.9 

1 

. 9982059 

.9939448  .9910294 

. 9887367 

200.0 

422.  1 

420.  B 

416.7  411.7 

409.  1 

.9466248 

. 9437094 

.9345145  .9233012 

.9174703 

250.0 

411.9 

409 

399.7  388.4 

352.6 

.9237497 

.9172460 

.8963893  .3710473 

. 8580399 

300.0 

410.5 

406.7 

393.7  375.5 

366.9 

.9206100 

.9120379 

.8829334  .8421171 

. 8228302 

Case  No 

Laminate 
Temperature 
(Deg.  F> !  0.00  ! 

.  241-260  C45., -45.328 

Nondi mensi onal  Time 

o.ooi  :  o.oi  !  o. i  : 

0.5 

80.0 

428.9 

428.7 

428.3 

427.7 

427.3 

1 

.9995337  . 

9986011 

.9972021 

. 9962695 

200.0 

408.4 

407 

402 

393.5 

389.5 

. 9522033 

.9439391  . 

9372814 

.9174633. 

.9081371 

250.0 

387 . 4 

384.  2 

372.  1 

351.1 

341.1 

.9032408 

.  8957799  . 

8675682 

.8186057 

. 7952903 

300.0 

381.7 

376.5 

q 

v»*  te/  w  ■  w 

319 

300.9 

.8899510 

.87/8270  . 

8295640 

. / 437631 

.7015621 
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Case  No.  401-420  (0.  , 45. , -45. , 90. > S 


Laminate  Nondimensional  Time 

Temperature 


(Deg.  F) 

:  o.oo  ' 

0.001  ! 

0.01  ! 

0.  1  ! 

0.5 

80.0 

514.8 

513.8 

511 

507.8 

506.8 

1 

. 9980575 

.9926185 

.9864025 

. 9844600 

200.0 

493.3 

490.8 

482.9 

474.5 

471.4 

. 9582362 

.9533800 

.9380342 

.9217172 

.9156954 

250.0 

482.3 

476.6 

458.7 

439.8 

432.6 

.9368687 

.9257964 

. G910256 

.8543124 

. 8403263 

250.0 

480.  4 

472.3 

445.9 

417.  1 

405 

.9331779 

.9174437 

.8661616 

.8102176 

.7867133 

Case 

No. 

421-440 

(90. ,45. ,-45. ,0. >S 

Laminate 

Temperature 

Nondi  i 

mensional  Time 

(Deg.  F)  !  0. 

00 

!  0.001 

:  o.oi  :  o.  i  ! 

0.5 

80.0 

445.9 

444.3 

440.3 

437.3 

436.3 

1 

.99641 18 

.9374411 

.9807132 

.9784705 

200.0 

422.  1 

419.4 

411.4 

403 

399.7 

.9466248 

. 9405696 

.9226284 

.9037901 

.8963893 

250.0 

411.9 

406 

387.8 

369 

361.2 

.9237497 

.9105131 

.8697017 

. 8275398 

.8100471 

250.0 

410.5 

403 

378 

349.6 

336.6 

.9206100 

.9037901 

.9477237 

. 7840323 

.7548778 

Case  No. 

Laminate 

Temperature 
(Deg.  F) !  0.00  ! 

441-460  (45. , -45. ) 2S 

Nondimensional  Time 

0.001  S  0.01  !  0.1 

0.5 

80.0 

428.9 

428.  4 

427.2 

425.8 

425.4 

1  . 

9908342  . 

9960364 

. 9927722 

.9719396 

200.0 

408.4 

406.5 

399.7 

386.8 

379.5 

.9522033  . 

9477734  . 

9319189 

.9018419 

.8848216 

250.0 

387,4 

383.  1 

367.3 

335.2 

316 

.9032403  . 

0932152  . 

8363760 

.7815342 

.7367605 

250.0 

33)  .7 

375.  1 

3*9.  1 

290.7 

245.7 

.8899510  . 

8745628  . 

8139426 

.6777804 

. 5728603 
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Case  No.  601-620  <0. , 45. , -A5. , 90. > S 
Radius  =  24  in. 


Laminate  Nandi  men si  on a 1  Time 

Temperature 


(Deg.  F) 

!  0.00  ! 

0.001  ! 

0.01  ! 

0.  1  ! 

0.5 

80.0 

290.8 

290.5 

209.6 

288.  5 

288.  1 

1 

. 998965 

.995373 

.992091 

.990715 

200.0 

279.3 

270.  7 

276.  8 

274.5 

273.  4 

. 960454 

. 95339 1 

.951857 

.943948 

.940165 

250.0 

273.5 

272.  1 

26e 

262.9 

260.4 

.940509 

. 935695 

.921596 

. 904035 

.895461 

300.0 

272.5 

270.7 

265 

257,7 

254 

.937070 

. 930330 

.911279 

.886176 

.873453 

Case  No,  701-720  <0. , 45. , -45. , 90. ) S 
Radius  *  48  in. 

Laminate  Nondimensicnal  Time 

Temperature 


(Deg.  F) 

1  0.00  ! 

0.001  ! 

H 

o 

« 

o 

0.  1  ! 

0.5 

80.0 

170.3 

170.2 

169.7 

169.  1 

163.9 

1 

.999413 

. 996477 

.992954 

.991779 

200.0 

164.5 

164.2 

163.2 

162 

161.4 

.965942 

.964181 

. 958309 

.951262 

. 947739 

250.0 

161.5 

160.8 

158.6 

156 

154.7 

.948326 

.944216 

.931293 

.916031 

. 90S397 

300.0 

160.9 

160 

157.  1 

153.3 

151,4 

.944803 

.939518 

.922490 

.900176 

.889019 

105 


Figure  38-  Case  No.  1-20;  (0. ,+45. ,-&5« *90. ) 
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Dimensionless  Time  (KT/n  ) 


Figure  42,  Case  No,  121-140;  (90. ,+45. , -45. ,0. ) 


Dimensionless  Time  (KT/h") 


Fig'sre  Case  No.  201-220;  (Q,  ,+45» « 


5., 90.) 


Dimensionless  Time  (KT/n  ) 


Figure  45.  Case  No.  221-240;  (90. ,+45* »-45« »0. ) 


m 


[#1 


m 


0.80  H 


0.75  H 


•  80  F 

•  200  F 
A  250  F 

m  300  f 


0.70  -j- 
0 


0.001  0.01  0.1 

Dimensionless  Time  (KT/h^) 


0 


Figure  $1.  Case  No.  701-720;  (0.  ,-t-45»  »-^5*  »9 0.; 
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CASE  NO .  40 

PLOT  NO.  6,  UNIT  I,  STEP  0 
DISPLACEMENT  CONTOURS 
RU  COMPONENT 

MODEL  SCRL.E  -  .5CCCE-GG.  CR!£\'7.  - 

SOL U T i C N  STRLE  =  ,6543E^07 


O.CC. 


G.  CC. 


CASE  NO .  4Q 

PLOT  NO.  7 ,  UNIT  I.  STEP  0 
DISPLACEMENT  CONTOURS 
RV  COMPONENT 

MODEL  SCPLE  --  .SDOGE-CG.  ORIENT.  -  C.CG.  C.CO. 
SOLUTION  SCP.LE  "  . iS^6E^C7 


CASE  NO .  60 

PLOT  NO .  1 ,  UNIT  U 

DISPLACEMENT  GEGMETR 
UVW  COMPONENT 


OOEL  SCALE  - 
OLUTICN  SCALE 


.SOOGE+OG.  C 
.  1  GGGE  +  G  4 


CASE  NO.  60 

PLOT  NO.  2,  UNIT  0,  STEP  0 
EIGENVECTOR  GEOMETRY 
UVW  COMPONENT  ,  MODE  1 

MODEL  SCALE  -  .5G00E-00.  CRIENT.  =  O.OD. 

SOLUTION  SCALE  -  .8000E-GC 


Figure  53.  Case  60,  UVW  Eigenvector  Component 


CASE  NO.  60 

PLOT  NO.  4,  UNIT  1.  STEP  0 
DISPLACEMENT  CONTOURS 
V  COMPONENT 

.  50CCE  +  00  .  0  "  i  E.  N  7  .  - 

.  5935E-G7 


MODEL  SCALE  = 
SOLUTION  SCALE 


0.00. 


0.00. 


CASE  NO.  60 

PLOT  NO.  5,  UNIT  1 ,  STEP  0 
DISPLACEMENT  CONTOURS 
W  COMPONENT 

MODEL  SCALE  =  .5C0GE+C0.  G=?!ENT.  -  G.OC.  O.GO. 

SOLUTION  SCALE  =  . 391  !EtQ6 


CASE  NO.  60 

PLOT  NO .  6,  UNIT  1,  5 

DISPLACEMENT  CONTOURS 
RU  COMPONENT 

MODEL  SCALE  -  .  SGOQE  +  GC.  GR  I  ENT 

SOLUTION  SCALE  -  .S9S3E<-L'6 


0 


CASE  NO.  60 

PLOT  NO,  1,  UNIT  1,  STEP 
DISPLACEMENT  CONTOURS 
RV  COMPONENT 

MGOEL  SCSLE  -  .EOGOE+OO.  ORIENT.  O.GO.  G.GO. 

SGLUTICN  SCGLE  =  ,SG79E<-G6 


Figure  63.  Case  60,  RV  Displacement  Component 
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*e  ro 


Case  120 #  V  Displacement  Conpcnen 


CASE  NO .  12G 

PLOT  NO.  b,  UNIT  i.  S 
DISPLACEMENT  CONTOURS 
W  COMPONENT 

MODEL  SCALE  =  . 5CCCL+00 .  CRIENI 

SOLUTION  SCALE  -  ,S348E‘0? 


CASE  N 0 .  120 
PLOT  NO.  6.  UN 
DISPLACEMENT  CO 
RU  COMPONENT 


MODEL.  SCALE  -  . 50G0E 

SOLUTION  SCALE  -  .42 


CL  Displacement  Component 


CASE  NO .  120 

PLOT  NO.  7 ,  UNIT  1,  S^EP  G 
DISPLACEMENT  CONTOURS 
RV  COMPONENT 

MODEL  SCALE  -  .5GCCE-G0.  0  R  I E N T .  =  0.00.  0.00, 

SOLUTION  SCALE  =  .SC2SE-C7 


Figure  69.  Case  120,  RY  Displacement  Component 
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Figure  ?2. 


CASE  NO .  140 

PLOT  NO.  7 .  UNIT  1 .  STEP  0 
DISPLACEMENT  CONTOURS 
RV  COMPONENT 

MODEL.  SCALE  -  .hCGOE^GO.  GR  !  EN  T .  -  0.00.  C.CC 

SOl.U T  I G'l  SCALE  -  .  M2CE*07 


CASE  NO .  16G 

PLOT  NO .  6.  UN[T  : ,  STEP  0 

DISPLACEMENT  CONTOURS 
RU  C  0  N  P  O  N  E N  T 

MODEL  SCRLE  -  . 5CCCE-0G.  ORIENT.  -  O.CC.  O.CO. 

SOLUTION  SC. RLE.  =  .  4672E«-06 


CASE  NO .  160 

PLOT  NO.  7,  UNIT  1,  STEP  0 
DISPLACEMENT  CONTOURS 
RV  COMPONENT 

MODEL  SCRLE  =  .SGGGE*00.  ORIENT.  =  O.CC 
SOLUTION  SCALE  -  .520^06 


Figure  81.  Cace  ioQ,  RV  Displacement  Ccr.por.ont 


Appendix  D 

Comparisons  of  Calculated  and  Predicted  Bifurcation  Loads 
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Case 

STAGS-C1 

Predicted 

No. 

N  Value 

N  Value 

X 

X 

21 

446.0 

446.00 

22 

445. 1 

445.69 

23 

442.9 

444.57 

24 

440.9 

441.48 

25 

440.  1 

439.79 

26 

422.2 

420.98 

27 

420.8 

420.34 

28 

417.0 

417.37 

29 

412.9 

411.15 

30 

410.6 

407.33 

31 

412.0 

410.24 

32 

409.0 

40e.G2 

33 

400.5 

403.43 

34 

392.  1 

388.67 

35 

386.7 

380.33 

36 

410.5 

408.68 

37 

406.8 

406.75 

38 

395.4 

399. 17 

39 

382.7 

378.51 

40 

375.3 

366.56 

41 

428.9 

428.90 

42 

428.6 

428.79 

43 

427.7 

428.39 

44 

426.7 

427.30 

45 

426.  4 

426.70 

46 

408.4 

404.65 

47 

407.8 

403.89 

48 

406.  1 

400.70 

49 

401.4 

392.06 

50 

397.7 

336.81 

51 

387.4 

379.29 

52 

386.3 

377.54 

53 

382.5 

370. 16 

54 

371.4 

350.23 

55 

362.2 

338.09 

56 

381.7 

372,58 

57 

380.  3 

369.74 

58 

375.0 

357, 63 

59 

356,6 

324.94 

60 

340.4 

304.86 

150 


Percent 

Error. 


.00 
.  13 
.38 
.  13 
-.07 
-.29 
-.  1  1 
.21 
-.42 
-.80 
-.43 
-.04 
.73 
-.87 
-1.65 
-.44 
-.01 
.95 
-1.09 
-2.33 


.00 
.04 
.  16 
.  14 
.07 
-.92 
-.96 
-1.33 
-2.33 
-2.74 
-2.09 
-2.27 
-3.22 
-5.70 
-6.66 
-2.39 
-2.78 
-4.63 
-8.88 
-10.44 


Case 

No. 


Predicted 
N  Value 
x 


Percent 

Error 


STAGS-C1 
N  Value 
x 


101 

502.80 

502.80 

.00 

102 

502.20 

502.54 

.07 

103 

500.50 

501.61 

.22 

104 

498.50 

499.03 

.  11 

105 

497.60 

497.62 

.00 

106 

481.20 

480.49 

-.  15 

107 

480.00 

479.90 

-.02 

108 

476.60 

477.63 

.22 

109 

472.40 

471.43 

-.20 

110 

470.30 

467.89 

-.51 

111 

470.30 

469.52 

-.  17 

112 

467.80 

468.22 

.09 

113 

460. 30 

463.21 

.63 

114 

451.00 

449.53 

-.33 

115 

'446.30 

441.71 

-1.03 

116 

468.50 

467.65 

-.  18 

117 

465.20 

465.83 

.  14 

118 

455.00 

458. 63 

.80 

119 

441.70 

439.03 

-.60 

120 

434.60 

427.58 

-1.61 

121 

444. 10 

444. lO 

.00 

122 

443.30 

443.79 

.  11 

123 

441.20 

442.60 

.  33 

124 

439.20 

439.61 

.09 

125 

438.40 

437.92 

-.11 

126 

420.60 

419. 16 

-.34 

127 

419.30 

418.51 

-.  19 

128 

415.50 

416.06 

.  13 

129 

411.50 

409 . 34 

-.52 

130 

409.20 

405.54 

-.90 

131 

410.50 

408.40 

-.51 

132 

407 . 70 

406.99 

-.  18 

133 

399.00 

401.60 

.65 

134 

389.90 

386.87 

-.78 

135 

385.00 

378.54 

-1.68 

136 

409.00 

406.83 

-.53 

137 

405.50 

404.90 

-.  15 

138 

393.80 

397.33 

.90 

139 

380.90 

376.69 

-1.11 

140 

373 . 50 

364.74 

-2.34 

Case 

STAGS-C1 

Predicted 

Percent 

No. 

N  Value 

N  Value 

Error 

X 

X 

141 

421.6 

421.60 

.00 

142 

421.3 

421.50 

.05 

143 

420.5 

421.11 

.  15 

144 

419.6 

420.06 

.11 

145 

419.2 

419.48 

.07 

146 

401.0 

397.32 

-.92 

147 

400.5 

396.55 

-.99 

148 

398.6 

393.34 

-1.32 

149 

393.6 

384.65 

-2.27 

150 

389.8 

379.36 

-2. 68 

151 

379.7 

371.62 

-2.  13 

152 

373.5 

369.86 

-2.28 

153 

374.3 

362.43 

-3.  17 

154 

362.3 

342.35 

-5.51 

155 

352.9 

330. 11 

-6.46 

156 

373.9 

364.80 

-2.44 

157 

372.3 

361.93 

-2.7? 

158 

366.0 

349.70 

-4.45 

159 

345.8 

316.71 

-8.41 

160 

328.9 

296.44 

-9.87 

401 

514.80 

514.80 

.00 

402 

513.80 

514.23 

.09 

403 

511.00 

512.39 

.27 

404 

507.80 

507.72 

-.02 

405 

506.80 

506.33 

-.09 

406 

493.30 

492.23 

-.22 

407 

490.80 

491.04 

.05 

408 

482.90 

486.45 

.73 

409 

474.50 

474.78 

.06 

410 

471.40 

470.05 

-.29 

411 

432.30 

481.10 

-.25 

412 

476.60 

478.47 

.39 

413 

458.70 

468.32 

2. 10 

414 

439 . 80 

442.53 

.  63 

415 

432.60 

432.28 

-.07 

416 

480.40 

479. 19 

-.25 

417 

472.30 

475.51 

.68 

418 

445.90 

460.91 

3.37 

419 

417.  10 

423.34 

1.50 

420 

405.00 

406. 15 

.28 

152 


STASS-C1 
N  Value 
x 


Predicted 
N  Value 
x 


Percent 

Error 


290.80 

290.50 
289.60 

288.50 
288. 10 
279.30 

278.70 

276.80 

274.50 

273.40 

273.50 
272. 10 

26e.00 

262.90 

260.40 

272.50 

270.70 
265.00 

257.70 
254.00 


170.30 

170.20 

169.70 

169. 10 

168.90 

164.50 

164.20 

163.20 
162.00 

161.40 

161.50 
160.80 
158.60 
156.00 

154.70 

160.90 
160.00 

157. 10 

153.30 

151.40 


290.80 

290.65 
290. 12 

288.65 

287.84 
278.05 
277.72 
276.42 
272.87 

270.85 
271.76 
271.02 
268. 15 
260.33 

255.85 
270.69 

269.65 
265.52 
254.31 
247.75 


170.30 

170.21 

169.90 
169.04 
168.57 
162.84 
162.64 
161.88 
159.80 
158.61 
159. 15 
158.72 
157.04 
152.45 
149.83 
158.52 

157.91 
155.50 
148.93 
145.09 


.00 
.05 
.  18 
.05 
-.09 
-.45 
-.35 
-.  14 
-.59 
-.93 
-.63 
-.40 
.06 
-.98 
-1.75 
-.67 
-.39 
.20 
-1.32 
-2.46 


.00 
.01 
.  12 
-.03 
-.20 
-1.01 
-.95 
-.81 
-1.36 
-1.73 
-1.45 
-1.30 
-.99 
-2.27 
-3.  15 
-1.48 
-1.31 
-1.02 
-2.85 
-4. 17 
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